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Figure  2.1.  Velocity  models  used  in  the  synthetic  seismogram  calcula-  8 

tions.  The  solid  curve  is  the  earth-flattening  transformation  for  a  homo¬ 
geneous  upper  mantle.  Note  the  break  in  the  velocity  axis. 

Figure  2.2.  Temporal  Pn  displacement  amplitude  (1-15  Hz)  as  a  func-  9 

tion  of  epicentral  distance  for  the  three  elastic  earth  models.  The  solid 
curve  is  the  analytic  attenuation  of  a  canonical  head  wave  scaled  to  the 
calculated  amplitude  for  Model  1  at  200  km. 

Figure  2.3.  Record  section  plot  of  Pn  particle  velocity  for  (a)  Model  1,  11 

(b)  Model  2  and  (c)  Model  3.  Travel  times  are  reduced  at  8.1  km/s. 

The  amplitudes  of  each  trace  are  scaled  independently. 

Figure  2.4.  Perspective  plot  of  Pn  displacement  spectra  as  a  function  of  14 

epicentral  distance  computed  for  ( a )  Model  1,  {b)  Model  2  and  (c) 

Model  3. 

Figure  2.5.  ( a )  Schematic  ray  diagram  for  the  primary  turning  ray  ( P )  17 

and  the  first  whispering  gallery  phase  (WG1).  (b)  Record  section  plot  of 
Pn  particle  velocity  for  Model  3  reducted  at  8.1  km/s.  Travel-time 
curves  for  P,  WG 1,  and  the  head  wave  are  superimposed. 

Figure  2.6.  Temporal  Pn  amplitude  (1-15  Hz)  for  Models  2  and  4.  The  19 

solid  curve  is  the  analytic  attenuation  of  a  canonical  head  wave. 

Figure  2.7.  Record  section  plot  of  Pn  particle  velocity  for  Model  4.  20 

Travel  times  are  reduced  at  8.1  km/s.  The  amplitudes  of  each  trace  are 
scaled  independently. 

Figure  2.8.  Perspective  plot  of  Pn  displacement  spectra  as  a  function  of  21 

epicentral  distance  computed  for  Model  4. 

Figure  3.1.  Map  showing  the  station  locations  (triangles)  and  event  30 

locations  (numbers)  for  the  NRDC  data  set.  The  numbers  correspond  to 
the  event  numbers  in  Table  3.2. 

Figure  3.2.  Pn  geometric  spreading.  The  solid  curve  is  the  Pn  attenua-  34 

tion  with  distance  at  1.25  Hz  calculated  for  an  elastic  earth  model  with 
upper  mantle  velocity  gradient  equivalent  to  the  earth-flattening  transfor¬ 
mation  of  a  homogeneous  medium  (Model  2,  Section  2).  The  dashed 
curve  plots  a  power  law  distance  dependence  r-1,25  determined  by  a 
least  squares  fit  to  the  synthetic  Pn  attenuation. 
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Figure  3.3.  Lg  Q(f)  models  consistent  with  the  NRDC  observations. 

The  value  of  Lg  Q  at  1  Hz  ( Q0 )  is  plotted  against  the  exponent  of  a 
power  law  frequency  dependence  (n).  These  models  combine  with 
different  source  levels  to  produce  data  variances  that  differ  by  less  than 
1%. 

Figure  3.4.  Derived  source  levels  (S0)  for  (a)  Lg  Q(f)  =  650  and  (b)  Lg  37 

Q(f)  =  345/ 04.  The  event  numbers  are  from  Table  3.2.  Events  recorded 
at  multiple  stations  have  an  S0  estimate  from  each  station.  Note  that  the 
higher  Q0  model  reduces  the  variance  in  the  derived  source  levels  from 
multiple  station  observations. 

Figure  3.5.  Moment  versus  magnitude  for  presumed  earthquakes  in  the  39 

NRDC  data  set.  The  moments  are  calculated  from  the  Lg  source  levels 
using  (3.5)  and  are  plotted  against  the  mb  reported  in  the  PDE  Bulletin 
(Events  1,  2,  3,  5,  and  9  in  Table  3.2).  The  solid  lines  are  moment 
versus  magnitude  relations  determined  in  previous  studies  (Table  3.3). 

The  dashed  horizontal  lines  are  moments  of  earthquakes  not  reported  in 
the  PDE  Bulletin. 

Figure  3.6.  Pn  Q(f)  models  consistent  with  the  NRDC  observations.  42 

The  value  of  Pn  Q  at  1  Hz  (Q0)  is  plotted  against  the  exponent  of  a 
power  law  frequency  dependence  (q).  These  models  combine  with 
different  source  levels  to  produce  data  variances  that  differ  by  less  than 
1%. 

Figure  3.7.  Derived  source  levels  (S0)  for  (a)  Pn  Q(f)  =  1175  and  (b)  Pn  43 

Q{f)  =  30Qfas.  The  event  numbers  are  from  Table  3.2.  Events  recorded 
at  multiple  stations  have  an  S0  estimate  from  each  station. 

Figure  3.8.  Moment  versus  magnitude  for  presumed  explosions  in  the  44 

NRDC  data  set.  The  moments  arc  calculated  from  the  Pn  source  levels 
using  (3.6)  and  are  plotted  against  Lg  magnitude  (Events  11,  13,  14,  16, 

18,  19-21  in  Table  3.2).  The  asterisks  are  moments  calculated  for  the 
Pn  Q(J)  =  3 QOf05  model  and  the  squares  arc  for  the  Pn  Q(f)  =  1175 
model.  Note  that  the  moments  for  the  lower  Q0  model  are  about  a  fac¬ 
tor  of  2  higher  than  the  moments  associated  the  higher  Q0  model.  The 
solid  line  is  the  explosion  moment  versus  magnitude  relation  determined 
by  Sereno  et  al.  [1988]  from  NORESS  observations. 

Figure  3.9.  Ambient  noise  displacement  spectra  for  a  5-s  window.  The  46 

solid  curve  is  the  average  noise  spectrum  at  the  three  NRDC  stations 
calculated  from  26  samples  taken  prior  to  Pn.  The  dashed  curves  indi¬ 
cate  one  standard  deviation.  The  dotted  curve  is  the  average  ambient 
noise  spectrum  at  NORESS  based  on  43  samples  taken  prior  to  Pn 
[Henson  and  Bache,  1988]. 
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Figure  3.10.  Theoretical  Lg  spectra  at  800  km  for  a  magnitude  3.0 
event  based  on  the  inversion  results  using  the  NRDC  (solid)  and 
NORESS  data  (dashed).  The  bottom  curve  is  the  average  NORESS 
ambient  noise  scaled  to  the  Lg  signal  window  length  at  800  km.  This 
curve  is  plotted  only  for  reference.  The  actual  Lg  spectrum  is  contam¬ 
inated  by  the  coda  of  previous  arrivals  which  is  larger  than  the  ambient 
noise. 

Figure  3.11.  Theoretical  Pn  spectra  at  800  km  for  a  magnitude  3.0  49 

explosion  based  on  the  inversion  results  using  the  NRDC  (solid  curves) 

and  NORESS  data  (dashed).  The  results  for  two  NRDC  Pn  models  are 

plotted  since  parameter  trade-offs  cannot  be  resolved  on  the  basis  of 

source  moment.  Model  parameters  are  listed  in  Table  3.4.  The  bottom 

curve  is  the  average  NORESS  ambient  noise  scaled  to  a  10-s  window 

length. 

Figure  A.l.  Waveforms  for  Event  1  (Table  3.2).  This  is  an  mb  4.3  62 

earthquake  reported  in  the  PDE  Bulletin.  The  calibration  values  are 
0.012781  and  0.012929  nanometers  per  digital  count  for  stations  KKL 
and  BAY,  respectively. 

Figure  A. 2.  Waveform  for  Event  2  (Table  3.2).  This  is  an  mb  4.3  63 

earthquake  reported  in  the  PDE  Bulletin.  The  calibration  value  is 
0.012929  nanometers  per  digital  count  for  station  KKL. 

Figure  A.3.  Waveform  for  Event  3  (Table  3.2).  This  is  an  mb  4.6  64 

earthquake  reported  in  the  PDE  Bulletin.  The  calibration  value  is 
0.012929  nanometers  per  digital  count  for  station  KKL. 

Figure  A.4.  Waveform  for  Event  4  (Table  3.2).  The  event  type  is  unk-  65 

nown.  The  calibration  value  is  0.012929  nanometers  per  digital  count 
for  station  KSU. 

Figure  A. 5.  Waveforms  for  Event  5  (Table  3.2).  This  is  an  mb  4.6  66 

earthquake  reported  in  the  PDE  Bulletin.  The  calibration  values  are 
0.411528  and  1.638321  nanometers  per  digital  count  for  stations  KKL 
and  BAY,  respectively. 

Figure  A.6.  Waveforms  for  Event  6  (Table  3.2).  The  event  type  is  67 

unknown.  The  calibration  values  are  0.012781  and  0.012929  nanome¬ 
ters  per  digital  count  for  stations  KKL  and  BAY,  respectively. 

Figure  A.  7.  Waveform  for  Event  7  (Table  3.2).  This  is  probably  an  68 

earthquake  but  was  not  reported  in  the  PDE  Bulletin.  The  calibration 
value  is  0.012929  nanometers  per  digital  count  for  station  BAY. 
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Figure  A.8.  Waveform  for  Event  8  (Table  3.2).  The  event  type  is  unk¬ 
nown.  The  calibration  value  is  0.006406  nanometers  per  digital  count 
for  station  KKL. 

Figure  A.9.  Waveforms  for  Event  9  (Table  3.2).  This  is  an  mb  4.4 
earthquake  reported  in  the  PDE  Bulletin.  Note  the  exceptionally  strong 
Sn  at  KKL  and  BAY.  The  calibration  value  is  0.006406  nanometers  per 
digital  count  for  all  three  stations. 

Figure  A.  10.  Waveforms  for  Event  10  (Table  3.2).  This  is  probably  an 
an  earthquake  but  was  not  reported  in  the  PDE  Bulletin.  The  magnitude 
of  this  event  is  approximately  3.3.  The  data  recorded  at  Bayanaul  are 
plotted  in  two  panels  because  there  is  a  data  gap  near  20  s.  Note  the 
strong  Sn  phase  at  both  stations.  The  calibration  value  is  0.006406 
nanometers  per  digital  count  for  both  stations. 

Figure  A.ll.  Waveforms  for  Event  11  (Table  3.2).  This  is  probably  an 
explosion.  The  Lg  magnitude  is  approximately  2.6.  Note  the  small 
Lg/Pn  ratio  for  this  event.  The  calibration  value  is  0.006406  nanometers 
per  digital  count  for  both  stations. 

Figure  A.  12.  Waveform  for  Event  12  (Table  3.2).  The  event  type  is 
unknown.  The  calibration  value  is  0.006406  nanometers  per  digital 
count  for  KSU. 

Figure  A.  13.  Waveforms  for  Event  13  (Table  3.2).  This  is  an  explo¬ 
sion  in  the  Lake  Balkash  area  [Thurber  et  al„  1989].  Note  the  presence 
of  a  strong  Rg  phase  at  stations  KSU  and  BAY.  The  Lg  magnitude  of 
this  event  is  approximately  2.5.  The  calibration  value  is  0.006406 
nanometers  per  digital  count  at  all  three  stations. 

Figure  A.  14.  Waveforms  for  Event  14  (Table  3.2).  This  is  probably  an 
explosion.  The  Lg  magnitude  of  this  event  is  approximately  2.3.  The 
calibration  value  is  0.006406  nanometers  per  digital  count  at  both  sta¬ 
tions. 

Figure  A.  15.  Waveforms  for  Event  15  (Table  3.2).  The  event  type  is 
unknown.  The  calibration  values  are  1.638321  and  0.006406  nanome¬ 
ters  per  digital  count  for  stations  KSU  and  KKL,  respectively. 

Figure  A.  16.  Waveform  for  Event  16  (Table  3.2).  The  event  is  prob¬ 
ably  an  explosion.  The  Lg  magnitude  for  this  event  is  approximately 
2.6.  The  calibration  value  is  0.006406  nanometers  per  digital  count  for 
station  KSU. 

Figure  A.  17.  Waveforms  for  Event  17  (Table  3.2).  This  event  is  an 
earthquake  in  the  Tien  Shan  area  [Thurber  et  al .,  1989].  The  calibration 


value  is  0.006406  nanometers  per  digital  count  for  all  three  stations. 

Figure  A.  18.  Waveforms  for  Event  18  (Table  3.2).  This  event  is  prob-  79 

ably  an  explosion.  The  Lg  magnitude  is  approximately  2.5.  The  cali¬ 
bration  value  is  0.006406  nanometers  per  digital  count  for  both  stations. 

Figure  A. 19.  Waveforms  for  Event  19  (Table  3.2).  This  event  is  a  10-  80 

ton  calibration  explosion  [Thurber  et  al. ,  1989;  Given  et  al.,  1989].  The 
vertical  scale  is  adjusted  for  the  waveform  at  KSU  (the  data  at  BAY  and 
KKL  are  not  clipped).  The  Lg  magnitude  for  this  event  is  approxi¬ 
mately  2.9.  The  calibration  value  for  station  KSU  is  0.006406  nanome¬ 
ters  per  digital  count  and  it  is  1.638321  for  stations  BAY  and  KKL. 

Figure  A.20.  Waveforms  for  Event  20  (Table  3.2).  This  event  is  a  10-  81 

ton  calibration  explosion  [Thurber  et  al.,  1989;  Given  ei  al.,  1989].  The 
vertical  scale  is  adjusted  for  the  waveform  at  KSU  (the  data  at  BAY  and 
KKL  are  not  clipped).  The  Lg  magnitude  for  this  event  is  approxi¬ 
mately  2.9.  The  calibration  values  are  0.050629,  0.102188,  and 
1.638321  nanometers  per  digital  count  for  stations  BAY,  KKL  and 
KSU,  respectively. 

Figure  A. 21.  Waveforms  for  Event  21  (Table  3.2).  This  event  is  an  82 

explosion  in  the  Lake  Balkash  area  (compare  to  Event  13).  Note  the 
strong  Rg  phase.  The  Lg  magnitude  for  this  event  is  approximately  2.9. 

The  calibration  value  is  0.012929  nanometers  per  digital  count  for  sta¬ 
tion  BAY. 
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1.  INTRODUCTION 


1.1  Project  Objectives 

The  objective  of  this  two-year  study  is  to  simulate  detection  and  location  capability  of 
seismic  networks  including  regional  stations  and  arrays  in  and  around  the  Soviet 
Union.  Three  specific  research  tasks  are: 

(1)  Enhance  and  validate  the  extended  version  of  the  Seismic  Network  Assessment 
Program  for  Detection  [SNAP/D,  Ciervo  et  al.,  1985]  called  SNAP/DX  [Bratt  et 
al.,  1987]  to  more  accurately  represent  the  treaty  monitoring  capability  of  seismic 
networks  including  regional  stations  and  arrays. 

(2)  Normalize  SNAP/DX  to  the  observed  performance  of  existing  stations  and 
expected  conditions  in  and  around  the  Soviet  Union. 

(3)  Apply  the  normalized  simulation  methods  to  asses  the  treaty  monitoring  capability 
of  existing  and  proposed  networks. 

1.2  Current  Status 

Nearly  half  of  the  research  effort  during  the  first  six  months  of  this  contract  was 
directed  toward  the  first  task  ( SNAP/DX  enhancement).  This  includes  incorporating  the 
ability  to  simulate  detection  capability  normalized  by  frequency  dependent  amplitude- 
distance  curves,  source  functions,  and  station  noise.  However,  since  formal  test  and 
evaluation  of  the  new  simulation  program  has  not  been  completed,  I  will  defer  its 
description  until  the  first  annual  report.  The  remainder  of  the  research  effort  during 
this  period  was  applied  to  the  normalization  task. 

T'  e  key  to  meaningful  simulations  of  detection  and  location  capability  is  accurate  nor¬ 
malization  to  the  observed  performance  of  existing  stations  and  arrays.  This  normali¬ 
zation  includes  the  determination  of  frequency  dependent  amplitude -distance  curves  for 
regional  seismic  phases,  noise  spectra  appropriate  for  primary  and  secondary  phases, 
source  spectra  and  scaling  relations,  and  beam  gain  for  array  stations.  Two  separate 
databases  of  regional  seismograms  and  spectra  are  being  assembled  as  part  of  this  nor¬ 
malization  effort.  These  are: 

•  NORESS/ARCESS.  Data  recorded  by  the  NORESS  and  ARCESS  arrays  in 
Scandinavia  from  regional  events  with  ML  between  2.0  and  4.0  will  be  used  as  the 
primary  basis  for  the  normalization.  NORESS  and  ARCESS  are  prototype  arrays 
for  regional  treaty  monitoring  and  are  within  regional  distances  to  parts  of  the 
western  USSR.  Also,  many  of  the  events  are  small  chemical  explosions  located 
in  Scandinavia  and  in  the  western  USSR.  Therefore,  these  data  provide  a  reason¬ 
able  basis  for  normalizing  the  detection  capability  of  a  network  of  stations  and 
arrays  in  and  around  the  Soviet  Union. 
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•  NRDC.  Data  collected  by  three  stations  in  Eastern  Kazakhstan,  USSR,  as  part  of 
a  joint  project  between  the  Natural  Resources  Defense  Council  (NRDC)  and  the 
Academy  of  Sciences  of  the  Union  of  Soviet  Socialist  Republics  are  used  to  com¬ 
pare  regional  wave  propagation  characteristics  and  noise  spectra  at  one  site  in  the 
USSR  with  those  in  Scandinavia.  The  attenuation  models  used  to  normalize  our 
simulations  are  derived  primarily  from  data  recorded  at  NORESS  and  ARCESS. 
The  NRDC  data  provide  an  excellent  opportunity  to  test  the  validity  of  extrapolat¬ 
ing  these  models  to  wave  propagation  in  the  Soviet  Union. 

We  have  processed  data  from  over  forty  regional  events  recorded  by  the  NORESS  and 
ARCESS  arrays  and  have  acquired  waveform  data  from  an  additional  sixty  events. 
The  processing  involves  interactive  picking  of  regional  phases,  phase  velocity  and 
azimuth  determination,  amplitude,  dominant  frequency,  and  the  calculation  of  Fourier 
signal  and  noise  spectra.  These  data  will  supplement  the  database  of  nearly  200  events 
recorded  at  NORESS  used  by  Sereno  et  al.  [1988]  to  derive  frequency  dependent 
amplitude-distance  curves  for  Pn  and  Lg.  Using  data  from  multiple  stations  should 
help  resolve  the  trade-offs  among  source  and  path  parameters  described  by  Sereno  et 
al.  [1988], 

Data  from  21  regional  events  recorded  by  one  or  more  of  the  NRDC  stations  have 
been  processed.  The  Pn  and  Lg  spectra  are  inverted  for  source  and  attenuation  param¬ 
eters  and  the  results  are  compared  to  models  derived  previously  using  NORESS  data 
[Sereno  et  al.,  1988].  Although  the  amount  of  NRDC  data  is  insufficient  to  resolve 
parameter  trade-offs  with  much  certainty,  the  spectra  are  consistent  with  the  attenua¬ 
tion  models  derived  for  Scandinavia.  The  NRDC  data  and  the  results  of  this  analysis 
are  presented  in  a  section  of  this  report. 

As  is  the  case  for  many  attenuation  estimation  methods,  the  generalized  inverse  tech¬ 
nique  used  by  Sereno  et  al.  [1988]  to  parameterize  the  distance  dependence  of  regional 
wave  spectra  requires  the  assumption  of  the  geometric  spreading  function.  This 
assumption  for  Lg  is  based  on  an  interpretation  as  a  sum  of  higher-mode  surface  waves 
and  has  been  carefully  studied  using  synthetic  seismogram  techniques  [Herrmann  and 
Kijko\  1983]  and  by  comparing  moments  derived  from  surface  waves  to  long-period  Lg 
amplitudes  [Street  et  al.,  1975],  However,  much  less  work  has  been  done  on  the 
geometric  spreading  of  Pn.  Spectral  ratio  methods  assume  that  Pn  spreading  is 
independent  of  frequency.  Other  methods  assume  a  power-law  range  dependence.  In 
this  report  I  investigate  the  sensitivity  of  Pn  geometric  spreading  to  the  velocity  gra¬ 
dient  in  the  upper  mantle.  Particular  emphasis  is  placed  on  the  frequency  dependence 
of  Pn  geometric  spreading  and  implications  for  Q  estimation  methods  that  assume 
frequency-independent  spreading. 

1.3  Outline  of  the  Report 

This  report  is  divided  to  two  major  sections.  The  first  section  reports  on  (1)  the  sensi¬ 
tivity  of  Pn  geometric  spreading  to  the  velocity  gradient  in  the  upper  mantle  deter¬ 
mined  using  synthetic  seismogram  techniques  and  (2)  the  relationship  between 


apparent  Pn  Q  determined  using  spectral  ratios  and  upper  mantle  anelasticity.  The 
second  section  reports  on  regional  data  recorded  in  eastern  Kazakhstan  as  part  of  the 
NRDC  experiment.  The  attenuation  of  Lg  and  Pn  are  parameterized  using  generalized 
inverse  techniques  and  the  results  are  compared  to  our  earlier  results  for  data  recorded 
at  NORESS  [Sereno  et  al.,  1988]. 
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2.  NUMERICAL  MODELING  OF  Pn  GEOMETRIC  SPREADING 

2.1  Introduction 

The  attenuation  of  regional  Pn  phases  has  important  implications  for  the  velocity  and 
anelastic  structure  of  the  upper  mantle  and  for  many  practical  issues  in  nuclear  explo¬ 
sion  monitoring  such  as  event  detection  and  identification.  Yet,  very  few  studies  of 
regional  Pn  attenuation  have  been  conducted.  While  it  has  long  been  recognized  that 
the  attenuation  of  Pn  is  a  strong  function  of  the  velocity  gradient  in  the  upper  mantle 
[e.g..  Hill,  1972],  there  have  been  no  quantitative  descriptions  of  the  sensitivity  of  Pn 
attenuation  to  upper  mantle  elastic  and  anelastic  structure  for  distances  between  200 
and  1000  km  and  frequencies  between  1  and  15  Hz  (ranges  relevant  for  regional 
nuclear  explosion  monitoring).  The  objectives  of  this  study  are  to  (1)  determine  the 
sensitivity  of  Pn  geometric  spreading  (between  200  and  1000  km  and  1  to  15  Hz)  to 
the  compressional-wave  velocity  gradient  in  the  upper  mantle  and  (2)  determine  the 
relationship  between  upper  mantle  anelasticity  and  apparent  Pn  Q  estimates  that  are 
based  on  the  assumption  of  frequency-independent  geometric  spreading. 

2.2  Approach 

Synthetic  seismograms  are  computed  for  three  elastic  models  with  different  velocity 
gradients  in  the  upper  mantle.  The  synthetic  seismograms  are  computed  using 
wavenumber  integration  [ Apsel ,  1979].  Synthetic  Pn  amplitude  spectra  for  these 
different  models  are  compared  as  a  function  of  distance  to  determine  the  sensitivity  of 
Pn  geometric  spreading  to  the  velocity  gradient  in  the  upper  mantle. 

To  determine  the  relationship  between  upper  mantle  anelasticity  and  apparent  Pn  Q 
estimates,  synthetic  Pn  phases  are  computed  for  an  anelastic  model  of  the  crust  and 
upper  mantle.  The  Pn  attenuation  for  this  model  is  compared  to  the  Pn  attenuation  for 
the  elastic  model  with  the  same  velocity  structure.  Apparent  Pn  Q  is  estimated  from 
the  synthetics  using  a  spectral  ratio  method  and  is  compared  to  the  known  Q  structure 
of  the  model. 

2.3  Results  for  Elastic  Earth  Models 

This  section  reports  the  results  of  the  synthetic  seismogram  calculations  for  three  elas¬ 
tic  models  of  the  crust  and  upper  mantle  in  Scandinavia.  The  elastic  models  are 
approximated  by  setting  the  quality  factor  for  shear  waves  and  compressional  waves 
equal  to  50000  in  the  crust  and  upper  mantle. 

2.3.1  Crustal  and  Upper  Mantle  Velocity  Models 

The  crust  is  approximated  by  a  homogeneous  40-km-thick  layer  with  compressional- 
wave  velocity  equal  to  6.5  km/s,  shear-wave  velocity  equal  to  3.75  km/s,  and  density 
equal  to  2700  kg/m3.  The  three  models  for  the  upper  mantle  velocity  are  listed  in 
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Table  2.1  and  compared  in  Figure  2.1.  The  velocity  at  the  top  of  the  mantle  (aMOHo) 
is  8.1  km/s  for  all  three  models. 

Model  1  is  a  simple  crustal  layer  over  a  mantle  halfspace.  This  model  is  included  to 
verify  the  accuracy  of  the  numerical  method  since  the  attenuation  for  this  model  can 
be  calculated  analytically.  Models  2  and  3  approximate  upper  mantle  gradients  of 
0.0013  s-1  and  0.0020  s_1,  respectively.  The  gradient  for  Model  2  is  approximately 
equal  to  the  earth -flattening  transformation  of  a  homogeneous  upper  mantle  (the  solid 
curve  in  Figure  2.1).  The  gradient  for  Model  3  is  slightly  greater  than  the  gradient  in 
the  upper  60  km  of  the  mantle  reported  by  King  and  Calcagnile  [1976]  for  Fennoscan- 
dia  and  western  Russia  (0.0017  s"1).  Mantle  density  (p)  is  calculated  from  the 
compressional-wave  velocity  (a)  using  the  relation,  p  =  379  a  +  252. 

The  wavenumber  integration  method  requires  that  the  earth  model  consist  of  a  stack  of 
homogeneous  layers.  Listed  in  Table  2.1  are  the  number  of  layers  used  to  approximate 
the  linear  velocity  gradient,  the  thickness  of  the  mantle  layers,  and  the  total  thickness 
of  the  mantle  gradient  zone.  At  larger  depths  the  model  is  a  halfspace  with  velocity 
equal  to  the  velocity  at  the  base  of  the  gradient.  The  layer  thickness  required  to 
approximate  a  linear  velocity  gradient  depends  on  the  wavelength  at  the  ray  turning 
point  [e.g.,  Burdick  and  Orcutt,  1979;  Chapman  and  Orcutt,  1985].  All  three  models 
use  a  mantle  layer  thickness  of  0.8  km  which  is  approximately  1.5  times  the  Pn 
wavelength  at  15  Hz  (the  maximum  frequency  used  in  the  calculation).  The  total 
thickness  of  the  gradient  zone  was  selected  so  that  the  pre-critical  reflection  from  the 
bottom  of  the  stack  did  not  interfere  with  the  primary  turning  ray  at  the  largest  epicen- 
tral  distance  (e.g.,  to  avoid  model  truncation  phases).  That  is,  the  gradient  extends  to 
great  enough  depth  so  that  the  synthetic  seismograms  at  distances  less  than  or  equal  to 
1000  km  are  independent  of  the  structure  below  the  gradient. 

Synthetic  displacement  spectra  for  an  explosion  source  were  computed  between  0  and 
20  Hz  for  distances  between  200  and  1000  km.  Both  source  and  receiver  were  located 
at  the  free  surface.  The  response  was  calculated  over  a  phase  velocity  interval 
appropriate  for  Pn  (8.0  knVs  to  20.0  km/s).  A  cosine-squared  taper  was  applied 
between  15  and  20  Hz  to  the  synthetic  spectra  prior  to  transforming  to  the  time 
domain.  The  explosion  source  time  history  is  represented  by  a  pressure  step  function 
and  the  source  moment  was  set  to  4.0  x  1013  nt-m. 

2.32  Broadband  Pn  Amplitude 

Figure  2.2  plots  temporal  Pn  displacement  amplitude  as  a  function  of  distance  for  each 
of  the  three  elastic  models.  The  amplitude  variation  of  a  canonical  head  wave  with 
distance  is  r~xa  t^a  where  r  is  epicentral  distance  and  lh  is  the  length  of  the  head  wave 
segment  [e.g.,  Aki  and  Richards ,  1980].  The  agreement  between  this  analytic  expres¬ 
sion  (the  solid  curve  in  Figure  2.2)  and  the  synthetic  Pn  amplitudes  computed  for 
Model  1  (squares  in  Figure  2.2)  gives  confidence  in  the  accuracy  of  the  numerical 
method. 
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Table  2.1  Elastic  Upper  Mantle  Models 


Model 

Number 

®Mcho 

(km/s) 

Velocity 
Gradient  (s-1) 

Number  of 
Layers 

Layer 

Thickness  (km) 

Total 

Thickness  (km) 

1 

8.1 

0.0000 

1 

oo 

oo 

2 

8.1 

0.0013 

76 

0.800 

60.0 

3 

8.1 

0.0020 

76 

0.800 

60.0 

Depth  Below  Moho 


Figure  2.1.  Velocity  models  used  in  the  synthetic  seismogram  calculations.  The  solid 
curve  is  the  earth-flattening  transformation  for  a  homogeneous  upper  mande.  Note  the 
break  in  the  velocity  axis. 
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Amplitude 


Temporal  Pn  Amplitude 


Figure  2.2.  Temporal  Pn  displacement  amplitude  (1-15  Hz)  as  a  function  of  epicentral 
distance  for  the  three  elastic  earth  models.  The  solid  curve  is  the  analytic  attenuation 
of  a  canonical  head  wave  scaled  to  the  calculated  amplitude  for  Model  1  at  200  km. 


All  three  models  produce  similar  Pn  attenuation  between  200  and  300  km  (the  cross¬ 
over  distance  for  these  models  is  approximately  110  km).  However,  the  head  wave 
amplitude  at  1000  km  is  close  to  a  factor  of  20  smaller  than  the  turning  wave  ampli¬ 
tude  for  Model  2  and  almost  a  factor  of  45  smaller  than  turning  wave  amplitude  for 
Model  3.  Note  also  that  Pn  amplitude  increases  between  600  and  1000  km  for  both 
linear  gradient  models.  These  results  are  similar  to  those  of  Hill  [1972],  although  his 
calculations  were  restricted  to  distances  less  than  400  km  and  were  for  a  thinner  cru¬ 
stal  model.  It  is  clear  from  Figure  2.2  that  the  Pn  amplitude  is  extremely  sensitive  to 
the  velocity  gradient  in  the  upper  mantle. 

2.3.3  Frequency  Dependence  ofPn  Geometric  Spreading 

The  Pn  waveforms  computed  for  the  linear  gradient  models  are  similar  to  the  head 
wave  at  short  ranges  (within  a  few  hundred  kilometers  of  the  critical  distance)  but 
evolve  into  a  turning  wave  with  increasing  distance.  The  distance  over  which  this 
transition  occurs  depends  on  the  crustal  thickness  and  the  velocity  gradient  in  the  man¬ 
tle.  Since  the  turning  ray  has  a  time  dependence  that  is  the  derivative  of  the  head 
wave  [e.g.,  Aki  and  Richards,  1980],  the  change  in  the  Pn  waveform  with  distance 
implies  a  strong  frequency  dependence  in  the  geometric  spreading  for  upper  mantle 
models  characterized  by  linear  velocity  gradients. 

For  example.  Figures  2.3a-2.3c  display  synthetic  Pn  particle  velocity  for  the  three 
elastic  models  described  in  Section  2.3.1.  Velocity  is  plotted  since  the  head  wave  pro¬ 
duces  static  displacement.  Note  that  the  amplitudes  are  scaled  to  be  the  same  on  each 
trace.  The  synthetic  Pn  waveforms  are  delta  functions  at  all  ranges  for  Model  1  (Fig¬ 
ure  2.3a).  The  geometric  spreading  for  this  model  is  independent  of  frequency.  At 
distances  less  than  about  300  km,  the  Pn  travel  paths  for  the  two  gradient  models  are 
almost  identical  to  the  Pn  head  wave  (the  turning  depth  is  less  than  1.5  km  below  the 
Moho).  Thus,  at  these  distances  the  Pn  waveforms  computed  for  the  gradient  models 
are  similar  to  the  Pn  head  waves  computed  for  Model  1.  However,  the  turning  wave 
dominates  for  ranges  greater  than  500-600  km  and  the  Pn  waveform  has  a  time  depen¬ 
dence  that  is  the  derivative  of  the  head  wave. 

Figure  2.4a  -  2.4c  are  perspective  plots  of  the  Pn  displacement  spectra  as  a  function  of 
epicentral  distance  for  the  three  elastic  models.  The  spectra  are  calculated  for  a  1-s 
window  centered  on  Pn.  These  plots  combine  the  observations  made  from  Figures  2.2 
and  2.3  into  a  quantitative  description  of  the  frequency  and  range  dependence  of  Pn 
geometric  spreading.  Note  that  the  spectral  shape  (to-1)  is  constant  for  all  ranges  for 
Model  1  and  is  equal  to  the  source  pressure  spectrum.  The  Pn  spectra  for  the  two  gra¬ 
dient  models  fall  off  as  to-1  at  short  ranges,  but  are  nearly  flat  beyond  500-600  km  due 
to  the  increasing  contribution  of  the  turning  ray. 

The  scalloping  in  the  Pn  spectra  for  Model  3  beyond  700  km  (Figure  2.4c)  is  due  to 
the  separation  of  the  primary  turning  wave  and  the  first  whispering  gallery  phase 
[Menke  and  Richards,  1980].  Figure  2.5a  is  a  schematic  ray  diagram  illustrating  the 
ray  paths  for  the  primary  turning  ray  ( P )  and  the  first  whispering  gallery  phase  (WG1) 
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Figure  2.4.  Perspective  plot  of  Pn  displacement  spectra  as  a  function  of  epicentral  dis¬ 
tance  computed  for  (a)  Model  1,  ( b )  Model  2  and  (c)  Model  3. 
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Figure  2.4.  Continued. 
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Figure  2.5.  ( a )  Schematic  ray  diagram  for  the  primary  turning  ray  ( P )  and  the  first 
whispering  gallery  phase  (WG1).  (b)  Record  section  plot  of  Pn  particle  velocity  for 
Model  3  reducted  at  8.1  km/s.  Travel-time  curves  for  P,  WG\,  and  the  head  wave  are 
superimposed. 


which  includes  an  upward  reflection  at  the  Moho.  Figure  2.5 b  is  the  Pn  record  section 
for  Model  3  with  travel-time  curves  for  P,  WG1,  and  the  head  wave  superimposed. 
Note  the  appearance  of  WG\  beyond  about  700  km.  At  longer  ranges  the  time  separa¬ 
tion  between  P  and  WG\  increases  and  the  spectra)  modulation  moves  toward  lower 
frequencies. 

2.4  Results  for  an  Anelastic  Earth  Model 

This  section  reports  the  results  of  the  synthetic  seismogram  calculations  for  an  anelas¬ 
tic  model  of  the  crust  and  upper  mantle.  The  velocities  and  densities  are  the  same  as 
Model  2  described  in  the  previous  section.  The  P -  and  S-  wave  quality  factors  were 
set  to  1000  and  500  in  the  crust,  respectively  [Kennett,  1985].  The  mantle  quality  fac¬ 
tors  were  set  to  2000  and  1000  for  compressional  and  shear  waves,  respectively. 
These  values  are  simply  used  for  demonstration  purposes  and  are  not  necessarily 
representative  of  upper  mantle  anelasticity  in  Fennoscandia.  For  notational  conveni¬ 
ence  this  will  be  referred  to  as  Model  4. 

Figure  2.6  plots  temporal  Pn  amplitude  for  Models  2  and  4.  The  solid  curve  is  the 
analytic  attenuation  for  a  head  wave.  The  Pn  attenuation  for  the  anelastic  model  with 
a  linear  gradient  (Model  4)  is  almost  identical  to  the  Pn  attenuation  of  a  head  wave  in 
an  elastic  medium  for  distances  less  than  400  km.  Figure  2.7  is  a  record  section  of  Pn 
particle  velocity  and  Figure  2.8  is  the  perspective  plot  of  Pn  displacement  spectra  as  a 
function  of  distance  for  Model  4.  The  obvious  effect  of  including  anelastic  attenuation 
is  the  reduction  in  the  amount  of  high  frequencies  (compare  Figure  2.7  with  Figure 
2.3 b  and  Figure  2.8  with  Figure  2.46).  The  large  reduction  in  high  frequencies 
between  900  and  1000  km  is  due  to  scalloping  caused  by  interference  with  the  first 
whispering  gallery  phase. 

Evernden  et  al.  [1986]  report  that  Pn  attenuation  in  the  eastern  United  States  can  be 
described  by  r~ 2  geometric  spreading  and  a  constant  Q  of  9000.  Their  results  are 
based  primarily  on  amplitudes  measured  between  1  and  3  Hz.  Chun  et  al.  [1989] 
report  that  the  attenuation  of  Pn  waves  between  1  and  15  Hz  in  eastern  Canada 
increases  with  increasing  frequency.  Their  model  is  consistent  with  the  Evernden  et 
al.  [1986]  model  at  low  frequency  but  predicts  smaller  Pn  amplitudes  at  frequencies 
greater  than  a  few  Hertz.  Figure  2.8  shows  that  the  Pn  attenuation  for  Model  4  is  only 
weakly  dependent  on  frequency  (except  for  the  modulation  due  to  the  interference  of 
the  whispering  gallery  phase).  Therefore,  the  frequency  dependence  of  Pn  attenuation 
described  Chun  et  al.  [1989]  suggests  that  if  the  velocity  gradient  in  the  upper  mantle 
is  greater  than  zero,  then  the  Q  (crustal  or  upper  mantle)  in  eastern  North  America  is 
lower  than  the  value  used  in  Model  4. 

2.5  Implications  for  Q  Estimation 

It  is  obvious  from  Figures  2.3 b  and  2.3c  that  unless  the  velocity  gradient  in  the  upper 
50-100  km  of  the  mantle  is  negative,  the  Q  estimates  that  are  derived  using  spectral 
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Figure  2.6.  Temporal  Pn  amplitude  (1-15  Hz)  for  Models  2  and  4.  The  solid  curve  is 
the  analytic  attenuation  of  a  canonical  head  wave. 
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Figure  2.7.  Record  section  plot  of  Pn  particle  velocity  for  Model  4.  Travel  times  are 
reduced  at  8.1  km/s.  The  amplitudes  of  each  trace  are  scaled  independently. 
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Pn  Displacement  Spectra 
Model  4:  5a  /  8z  =  0.0013  S  '(Q) 


Figure  2.8.  Perspective  plot  of  Pn  displacement  spectra  as  a  function  of  epicentral 
tance  computed  for  Model  4. 
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ratio  methods  will  be  biased  high.  This  is  because  the  spectral  ratio  method  assumes 
that  the  change  in  the  Pn  spectrum  with  distance  is  due  to  anelasticity.  However,  the 
differentiation  of  the  Pn  waveform  several  hundred  kilometers  beyond  the  critical  dis¬ 
tance  enhances  the  high  frequencies  relative  to  low  frequencies.  Thus,  the  apparent  Q 
includes  both  the  geometric  effect  (which  increases  the  amount  of  high  frequencies) 
and  the  anelastic  effect  (which  reduces  the  amount  of  high  frequencies).  The  bias  in  Q 
is  larger  for  spectral  ratios  calculated  over  distances  including  the  transition  zone  from 
head  wave  to  turning  wave  (e.g.,  400  to  700  km)  than  for  ratios  calculated  over  a 
smaller  interval  outside  the  transition  region  (e.g.,  800  to  1000  km). 

Spectral  ratios  were  computed  from  the  synthetic  Pn  phases  to  quantify  the  Q  bias. 
The  logarithm  of  the  ratio  of  the  displacement  spectra,  A(f,r),  at  two  ranges  is 
expressed  as 


log 


A(fri) 

Atfjri) 


log 


G('i) 

G(ra) 


7i  loge  Pi  -  t£f 
Q 


(2.1) 


where  G(r)  is  geometric  spreading,  t  is  travel  time,  and  /  is  frequency.  If  geometric 
spreading  and  Q  are  assumed  to  be  independent  of  frequency,  then  (2.1)  is  the  equation 
of  a  straight  line  with  slope  proportional  to  CT1  and  intercept  equal  to  the  logarithm  of 
the  geometric  spreading  ratio  at  the  two  ranges.  If  we  further  assume  that  the 
geometric  spreading  has  a  power  law  range  dependence,  r~n,  then  the  intercept  of  the 
spectral  ratio  is  given  by 


(2.2) 

Table  2.2 a  lists  Q  values  and  geometric  spreading  exponents,  -n,  derived  from  spectral 
ratios  of  the  synthetic  Pn  phases  computed  for  Model  1.  Table  2.2 b  lists  the  analytic 
values  for  these  parameters.  In  this  case  geometric  spreading  is  independent  of  fre¬ 
quency  and  the  spectral  ratio  method  accurately  recovered  the  Pn  spreading  rate  and 
very  high  Q  values.  Table  2.3a  lists  Q  values  and  geometric  spreading  exponents 
derived  from  spectral  ratios  of  the  synthetic  Pn  phases  computed  for  Model  2.  Nega¬ 
tive  Q  values  are  obtained  because  of  the  increasing  contribution  of  the  turning  ray  and 
the  differentiation  of  the  Pn  waveform.  The  spreading  rates  are  similar  to  those  of  a 
head  wave  between  200  and  400  km  but  increase  rapidly  at  larger  distances  and 
become  positive  between  700  and  1000  km.  Table  2.3 b  lists  Q  values  and  geometric 
spreading  exponents  for  Model  4.  In  general  the  estimated  Q  values  are  negative.  The 
lower  Q  values  associated  with  the  spectrum  at  1000  km  are  due  to  the  spectral  hole  at 
that  distance  near  15  Hz  which  is  caused  by  interference  with  the  first  whispering  gal¬ 
lery  phase.  Aside  from  this  anomaly,  the  effect  of  assuming  a  frequency-independent 
Pn  geometric  spreading  function  is  to  greatly  overestimate  upper  mantle  Q  using  the 
spectral  ratio  method. 


log 
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200 

300 

400 

500 

600 

700 

800 

900 

1000 

200 

9217 

16052 

18231 

21250 

23518 

25299 

17601 

24660 

300 

-3.13 

62097 

35677 

37623 

38421 

38860 

20746 

32420 

400 

-2.95 

-2.69 

25029 

31430 

34089 

35535 

18308 

30028 

500 

-2.83 

-2.58 

-2.45 

42231 

41624 

41317 

17156 

31278 

600 

-2.75 

-2.52 

-2.40 

-2.35 

41033 

40875 

14322 

29373 

700 

-2.69 

-2.48 

-2.38 

-2.33 

-2.30 

40718 

10805 

26831 

800 

-2.65 

-2.45 

-2.35 

-2.31 

-2.28 

-2.26 

6229 

22923 

900 

-2.59 

-2.39 

-2.29 

-2.23 

-2.17 

-2.10 

-1.91 

-13644 

1000 

-2.58 

-2.40 

-2.30 

-2.26 

-2.22 

-2.19 

-2.15 

-2.41 

Table  2.2a.  Q  estimates  and  the  exponent  of  a  power  law  geometric  spreading  func¬ 
tion  determined  using  spectral  ratios  for  synthetic  Pn  phases  computed  for  Model  1 
(layer  over  a  halfspace).  The  Q  estimates  are  in  the  upper  right  diagonal  and  the 
spreading  exponents  are  in  the  lower  left  diagonal.  Entries  are  given  for  each  pair  of 
epicentral  distances.  For  example,  the  Q  estimated  from  the  spectral  ratio  between  500 
and  900  km  is  17156  and  the  spreading  exponent  is  -2.23. 


200 

300 

400 

500 

600 

700 

800 

900 

1000 

50000 

50000 

50000 

50000 

50000 

50000 

50000 

50000 

-3.21 

50000 

50000 

50000 

50000 

50000 

50000 

50000 

400 

-2.99 

-2.68 

50000 

50000 

50000 

50000 

50000 

50000 

500 

-2.87 

-2.59 

-2.48 

50000 

50000 

50000 

50000 

50000 

600 

-2.78 

-2.53 

-2.43 

-2.37 

50000 

50000 

50000 

50000 

700 

-2.72 

-2.49 

-2.39 

-2.34 

-2.30 

50000 

50000 

50000 

800 

-2.68 

-2.46 

-2.37 

-2.31 

-2.28 

-2.25 

50000 

50000 

900 

-2.64 

-2.43 

-2.34 

-2.29 

-2.26 

-2.24 

-2.22 

50000 

1000 

-2.61 

-2.41 

-2.33 

-2.28 

-2.25 

-2.22 

-2.21 

-2.19 

Table  2.2 b.  Q  values  and  exponents  of  a  power  law  geometric  spreading  function 
determined  analytically  for  a  head  wave  propagating  in  a  medium  with  Q  =  50000 
(compare  to  Table  2.2a). 
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200 

300 

400 

500 

600 

700 

800 

900 

1000 

-1605 

-1567 

-1469 

-1399 

-1369 

-1390 

-1502 

-1826 

-3.17 

-1531 

-1410 

-1342 

-1320 

-1354 

-1487 

-1863 

400 

-2.82 

-2.32 

-1306 

-1263 

-1262 

-1316 

-1478 

-1932 

500 

-2.59 

-2.13 

-1.89 

-1223 

-1241 

-1319 

-1528 

-2137 

600 

-2.40 

-1.96 

-1.70 

-1.46 

-1260 

-1373 

-1667 

-2628 

700 

-2.23 

-1.78 

-1.50 

-1.25 

-0.99 

-1508 

-1988 

-4117 

800 

-2.05 

-1.59 

-1.28 

-1.00 

-0.70 

-0.36 

-2915 

-30453 

900 

-1.84 

-1.35 

-1.01 

-0.67 

-0.32 

0.10 

0.62 

3606 

1000 

-1.58 

-1.04 

-0.64 

-0.23 

0.20 

0.72 

1.37 

2.21 

Table  2.3a.  Q  estimates  and  the  exponent  of  a  power  law  geometric  spreading  func¬ 
tion  determined  using  spectral  ratios  for  synthetic  Pn  phases  computed  for  Model  2. 
Entries  are  in  the  same  format  as  those  in  Table  2.2a. 


200 

300 

400 

500 

600 

700 

800 

900 

1000 

8603 

-80941 

-10492 

-6181 

-5648 

-6146 

-9675 

36524 

-2.80 

-7095 

-4973 

-3930 

-3994 

-4577 

-7145 

68096 

400 

-2.53 

-2.15 

-3828 

-3213 

-3486 

-4204 

-7155 

24616 

500 

-2.33 

-1.95 

-1.68 

-2768 

-3337 

-4346 

-9142 

9902 

600 

-2.16 

-1.78 

-1.52 

-1.31 

-4200 

-6076 

-39297 

4618 

700 

-1.96 

-1.56 

-1.25 

-0.97 

-0.56 

-10977 

12368 

2717 

800 

-1.75 

-1.32 

-0.97 

-0.64 

-0.21 

0.20 

3956 

1673 

900 

-1.51 

-1.03 

-0.63 

-0.23 

0.25 

0.75 

1.36 

1061 

1000 

-1.19 

-0.65 

-0.17 

0.31 

0.89 

1.52 

2.31 

3.37 

Table  2.3 b.  Q  estimates  and  the  exponent  of  a  power  law  geometric  spreading  func¬ 
tion  determined  using  spectral  ratios  for  synthetic  Pn  phases  computed  for  Model  4. 
Entries  are  in  the  same  format  as  those  in  Table  2.2a. 
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2.6  Discussion 


The  amplitude  of  Pn  is  extremely  sensitive  to  the  velocity  gradient  in  the  upper  mantle. 
For  example,  the  synthetic  Pn  amplitude  at  1000  km  is  20  times  greater  for  a  model 
with  a  weak  upper  mantle  gradient  (0.0013  s-1)  than  the  canonical  head  wave  ampli¬ 
tude  for  a  layer  over  a  halfspace.  This  gradient  corresponds  to  the  earth-flattening 
transformation  of  a  homogeneous  upper  mantle.  Thus,  earth  sphericity  is  sufficient  to 
cause  significant  departure  of  Pn  attenuation  from  that  of  a  canonical  head  wave. 

The  geometric  spreading  of  a  canonical  head  wave  is  independent  of  frequency.  The 
geometric  spreading  of  Pn  for  a  model  with  increasing  velocity  with  depth  is  strongly 
dependent  on  frequency.  The  Pn  waveform  for  a  model  with  a  weak  velocity  gradient 
(like  that  introduced  by  sphericity  alone)  is  similar  to  a  head  wave  at  distances  within 
a  few  hundred  kilometers  of  the  critical  distance,  but  has  the  time  dependence  of  a 
turning  wave  at  longer  ranges.  Thus,  the  Pn  waveform  at  long  ranges  is  the  derivative 
of  the  Pn  waveform  near  the  critical  distance.  The  distance  over  which  the  transition 
occurs  depends  on  the  crustal  thickness  and  the  upper  mantle  velocity  gradient.  For 
the  models  considered  here,  the  transition  occurs  between  about  400  and  700  km. 
Burdick,  et  al.  [1988]  show  that  Pn  waveforms  from  explosions  at  the  Nevada  Test 
Site  are  consistent  with  a  turning  ray  representation  which  suggests  that  there  is  a  posi¬ 
tive  velocity  gradient  in  the  upper  mantle  in  that  region. 

Many  methods  for  estimating  anelastic  attenuation  assume  that  geometric  spreading  is 
independent  of  frequency.  This  can  lead  to  an  overestimate  of  Pn  Q  if  the  velocity  in 
the  upper  50-100  km  of  the  mantle  increases  with  depth.  For  example,  it  is  shown  in 
this  report  that  Q  cannot  be  recovered  from  synthetic  Pn  phases  computed  for  a  very 
simple  model  of  the  upper  mantle  (a  linear  gradient)  using  a  spectral  ratio  method. 
Therefore,  accurate  estimates  of  the  velocity  gradient  of  the  upper  mantle  are  required 
before  meaningful  estimates  of  upper  mantle  anelasticity  can  be  made  from  Pn  obser¬ 
vations. 
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3.  ATTENUATION  OF  Pn  AND  Lg  PHASES  RECORDED  IN  EASTERN 
KAZAKHSTAN 

3.1  Introduction 

Sereno  and  Bratt  [1989]  use  attenuation  models  and  ambient  noise  spectra  derived 
from  data  recorded  by  the  NORESS  array  in  Norway  to  estimate  the  detection  capabil¬ 
ity  of  a  network  of  seismic  arrays  in  and  around  the  Soviet  Union.  These  estimates 
have  a  large  uncertainty  since  it  is  not  know  how  well  NORESS  represents  stations 
inside  the  Soviet  Union.  Data  recently  obtained  from  three  stations  inside  the  Soviet 
Union  provide  an  excellent  opportunity  to  check  how  well  NORESS  represents  condi¬ 
tions  at  these  sites.  The  three  seismic  stations  were  operated  in  eastern  Kazakhstan, 
USSR,  as  part  of  a  joint  experiment  involving  the  National  Resources  Defense  Council 
(NRDC)  in  the  United  States  and  the  Academy  of  Sciences  of  the  Union  of  Soviet 
Socialist  Republics. 

Spectra  of  Pn  and  Lg  phases  from  21  regional  events  recorded  by  one  or  more  of  these 
NRDC  seismic  stations  are  inverted  for  source  and  attenuation  parameters  using  the 
method  of  Sereno  et  al.  [1988].  The  average  ambient  noise  spectrum  for  the  NRDC 
stations  is  computed  from  time  windows  preceding  Pn.  This  noise  spectrum  and  the 
parameterization  of  regional  wave  attenuation  are  compared  with  results  determined 
previously  for  Scandinavia  using  NORESS  data  [ Sereno  et  al.,  1988].  This  report 
summarizes  the  major  conclusions  from  the  comparison.  The  data  from  the  NRDC 
stations  are  too  few  to  normalize  estimates  of  the  detection  threshold  in  the  Soviet 
Union,  but  the  comparison  with  our  NORESS  results  allows  an  assessment  of  the  vali¬ 
dity  of  using  NORESS  performance  to  normalize  estimates  of  the  detection  threshold. 

3.2  Event  and  Station  Locations 

Three  seismic  stations  were  installed  in  1987  within  200  km  of  the  principal  under¬ 
ground  nuclear  explosion  test  site  in  eastern  Kazakhstan  [Berger  et  al.,  1988],  Each 
station  is  equipped  with  two  3-component  surface  seismometers  and  a  3-component 
borehole  seismometer.  The  station  locations  and  sensor  depths  of  the  borehole 
seismometers  are  listed  in  Table  3.1.  Berger  et  al.  [1988]  give  a  complete  description 
of  the  NRDC  sites  and  instrumentation. 

The  data  used  in  this  study  consist  of  Pn  and  Lg  phases  recorded  by  one  or  more  of 
the  three  NRDC  stations  from  21  regional  events  at  distances  between  200  and  1300 
km.  Event  magnitudes  are  approximately  2.3  to  4.6.  Table  3.2  lists  the  location,  ori¬ 
gin  time,  and  magnitude  of  each  event  used  in  this  study.  The  station  locations  and 
event  epicenters  are  plotted  on  a  map  in  Figure  3.1.  Events  1,  2,  3,  5  and  9  are  earth¬ 
quakes  listed  in  the  Preliminary  Determination  of  Epicenters  (PDE),  a  bulletin  pro¬ 
duced  by  the  United  States  Geological  Survey  (USGS).  The  location  of  Event  9  is 
from  Thurber  et  al.  [1989]  who  used  travel  times  and  azimuths  determined  from  data 
recorded  at  the  three  NRDC  stations  to  obtain  a  location  40  km  east  of  the  PDE 
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Table  3.1  NRDC  Station  Location 


Name 

Latitude 

Longitude 

Elevation  (m) 

Senor  Depth  (m) 
(Borehole  Sensor) 

Karkaralinsk  (KKL) 

49.333 

75.383 

1000 

66 

Karasu  (KSU) 

49.950 

81.083 

420 

101 

Bayanaul  (BAY) 

50.817 

75.550 

600 

99 
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Table  3.2  Event  Parameters 


Event 

Julian 

Day 

Origin 

Time 

Latitude 

Longitude 

Depth 

(km) 

mb 

Event 

Type 

1 

136 

18:21:13.7 

39.339 

73.812 

33.0 

4.3 

EQ 

2 

120 

6:54:57.0 

39.733 

74.594 

10.0 

4.3 

EQ 

3 

122 

21:38:33.2 

39.569 

74.541 

33.0 

4.6 

EQ 

4 

123 

8:52:08.0 

51.512 

89.874 

5 

146 

13:44:19.9 

42.923 

78.063 

20.0 

4.6 

EQ 

6 

164 

2:19:13.2 

47.300 

79.473 

7 

162 

18:30:04.0 

42.842 

80.968 

EQ 

8 

232 

8:57:32.5 

50.374 

68.615 

9 

234 

0:21:50.7 

44.129 

85.363 

58.0 

4.4 

EQ 

10 

239 

7:52:21.5 

44.148 

83.133 

EQ 

11 

239 

8:05:26.2 

49.982 

80.832 

EX 

12 

239 

8:51:21.2 

54.665 

86.488 

13 

239 

9:38:34.8 

46.900 

77.389 

EX 

14 

239 

11:04:40.5 

50.905 

83.464 

EX 

15 

242 

4:57:11.3 

47.808 

81.658 

16 

243 

7:57:42.9 

55.739 

81.080 

EX 

17 

244 

3:44:38.8 

43.808 

85.948 

EQ 

18 

244 

11:58:21.6 

54.316 

85.648 

EX 

19 

245 

7:00:00.0 

50.281 

72.172 

0.0 

EX 

20 

246 

7:00:00.0 

50.281 

72.172 

0.0 

EX 

21 

149 

9:28:17.0 

46.855 

78.057 

EX 
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location.  Event  17  is  a  smaller  earthquake  that  occurred  in  the  same  region  as  Event  9 
(Tien  Shan),  but  was  not  reported  in  the  PDE  [Thurber  et  al.,  1989).  Event  13  is  a 
mining  explosion  near  Lake  Balkash  [Thurber  et  al.,  1989).  Events  19  and  20  are 
10-ton  calibration  shots  detonated  west  of  the  city  of  Karaganda  [Given  et  al.,  1989). 
The  locations  of  the  remaining  events  in  Table  3.2  were  determined  by  seismologists 
at  the  Center  for  Seismic  Studies  (CSS)  based  on  travel  time  picks  from  the  NRDC 
data.  Event  21  is  close  to  the  mining  explosion  in  the  Lake  Balkash  area  reported  by 
Thurber  et  al.  [1989]  and  occurred  at  a  similar  time  of  day.  Also,  the  data  recorded  at 
Bayanaul  display  a  strong  Rg  phase  which  is  indicative  of  shallow  depth.  Thus,  Event 
21  is  probably  a  mining  explosion  in  the  Lake  Balkash  area.  Events  7  and  10  are 
probably  earthquakes.  Their  locations  are  between  the  PDE  locations  for  Events  5  and 
9.  The  remainder  of  the  events  are  of  unknown  origin,  though  it  is  likely  that  many 
are  near-surface  explosions. 

3.3  Data  Processing 

The  results  to  be  described  were  obtained  by  analysis  of  data  from  the  vertical  com¬ 
ponent  borehole  sensor.  All  available  short-period  data  were  processed  with  the  pro¬ 
cedures  described  below,  but  the  borehole  data  are  emphasized  because  the  surface 
seismometers  are  typically  characterized  by  higher  noise  levels  (particularly  at  KSU) 
[Berger  et  al.,  1988).  Also,  the  vertical  component  data  provide  direct  comparison  to 
out  previous  results  for  NORESS.  Appendix  A  displays  the  recorded  waveforms  for 
all  of  the  events  used  in  this  study. 

The  waveforms  for  each  event-station  pair  are  filtered  over  four  frequency  bands,  1-3, 
3-5,  4-8,  and  8-16  Hz.  The  broadband  and  filtered  traces  are  displayed  with  theoreti¬ 
cal  arrival  times  (based  on  the  locations  in  Table  3.2)  for  interactive  picking  of 
regional  phases.  Fourier  spectra  are  computed  for  each  phase.  A  symmetric  10% 
cosine  taper  is  applied  to  a  time  window  starting  0.3  s  before  the  picked  arrival  time. 
The  duration  of  the  analyzed  segment  depends  on  the  phase  identification.  For  Pn  the 
duration  was  set  to  10  s,  and  a  fixed  group  velocity  window  of  3.6  to  3.0  km/s  was 
used  to  deteimine  the  duration  for  the  Lg  spectral  estimate.  Noise  spectra  for  each 
phase  are  calculated  for  a  5-s  window  starting  5.5  s  before  the  picked  arrival  time. 
The  noise  spectra  are  normalized  to  the  signal  window  length  using 


where  N(f)  is  the  normalized  noise  spectrum,  Ts  is  the  window  length  used  to  calculate 
the  signal  spectrum,  Tn  is  the  window  length  used  to  calculate  the  noise  spectrum  (5 
s),  and  n(f)  is  the  noise  spectrum  computed  using  a  window  length  of  Tn. 

The  spectra  were  corrected  for  the  instrument  response  and  converted  to  units  of  dis¬ 
placement.  The  NRDC  borehole  instrument  response  is  flat  to  velocity  between  0.2 
and  100  Hz  [Berger  et  al.,  1988).  The  calibration  values  were  obtained  for  each 
waveform  from  the  Center  for  Seismic  Studies  (CSS)  in  units  of  nanometers  per  digital 
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count  at  1  s.  These  calibrations  were  verified  for  a  subset  of  the  data  by  H.  Given  at 
the  Institute  of  Geophysics  and  Planetary  Physics  [personal  communication].  The  dis¬ 
placement  spectrum  is  related  to  the  uncorrected  spectrum  by 


A(/)  = 


X(f) 

2  Kfl0 


(3.2) 


where  A(f)  is  the  displacement  spectrum  in  nm-s,  X(f)  is  the  spectrum  of  the 
uncorrected  data  in  count-s,  /  is  frequency  in  Hz,  and  /0  is  the  instrument  response  to 
velocity  in  counts/(nm/s).  The  calibration  values  reported  by  CSS  are  expressed  in 
terms  of  /0  as 

calib  =  (3.3) 

2  K  Iq 


where  Tc  is  the  calibration  period  (Tc  =  1  s).  Therefore,  the  signal  and  noise  spectra 
were  corrected  for  the  instrument  response  by  multiplying  by  the  calibration  value  and 
dividing  by  frequency. 


3.4  Attenuation  and  Source  Parameters  for  Regional  Phases 


The  method  of  Sere  no  et  al.  [1988]  is  used  to  estimate  attenuation  and  source  parame¬ 
ters  from  the  Pn  and  Lg  spectra  recorded  by  the  NRDC  stations.  The  instrument- 
corrected  amplitude  spectra  are  parameterized  as 


A(f,r)  =  S(f)  G(r,ro)  exp 


iZH 

Q(f) 


(3.4) 


where  S(f)  is  the  source  spectrum,  G(r,ro)  is  the  geometric  spreading  function,  and  the 
last  term  is  the  apparent  attenuation  for  travel  time  r.  The  spectra,  A(f,r),  are  corrected 
for  an  assumed  geometric  spreading  function  and  inverted  using  damped  least  squares 
for  two  parameters  describing  the  source  spectrum  and  two  parameters  characterizing  a 
power  law  frequency  dependence  of  Q(f). 


The  source  spectrum  is  assumed  to  falloff  as  f~2  beyond  a  comer  frequency  that  scales 
inversely  with  the  cube  root  of  the  long-period  level,  S0.  The  source  parameters 
estimated  by  the  inversion  are  the  long-period  source  level  for  each  event  (which  is 
proportional  to  seismic  moment)  and  a  single  parameter  that  relates  comer  frequency 
and  long-period  level  for  all  events.  The  relationship  between  source  moment  and  Lg 
S0  for  earthquakes  is 


SlHLg)  = 


4rcpcPc 


(3.5) 
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where  is  the  earthquake  moment,  pe  is  the  average  crustal  density,  and  (3e  is  the 
average  shear-wave  velocity  of  the  crust.  The  relationship  between  Pn  S0  and  earth¬ 
quake  moment  depends  on  the  source  radiation  pattern.  Since  the  focal  mechanisms 
for  the  events  studied  are  unknown,  it  is  not  possible  to  obtain  accurate  estimates  of 
source  moment  from  the  derived  Pn  source  levels.  However,  explosion  moments  can 
be  estimated  from  Pn  S0  using 


S$p(Pn)  = 


MS*? 


(3.6) 


where  ps  and  a,  are  the  near-surface  density  and  compressional-wave  velocity,  respec¬ 
tively.  Since  the  excitation  function  for  Lg  is  unknown  for  near-surface  explosions  it 
is  not  possible  to  estimate  explosion  moments  from  Lg  S0. 


In  this  report  the  geometric  spreading  function  is  assumed  to  be  independent  of  fre¬ 
quency  and  is  characterized  by  a  power  law  range  dependence 

G(r,r0)  =  (i /r)  for  r  <  r0 

G(r,r0 )  =  rjj1  {r^r)m  for  rtr0  ^7) 

where  r  is  epicentral  distance  and  r0  is  a  transition  distance  from  spherical  spreading  to 
spreading  rate  m.  The  Lg  phase  consists  of  higher-mode  surface  waves  which  are 
accurately  described  at  long  ranges  by  cylindrical  spreading  (m  =  1/2),  provided  the 
window  length  is  sufficient  to  encompass  the  entire  dispersed  wave  train.  By  compar¬ 
ing  the  long-period  amplitude  spectrum  of  Lg  to  moments  calculated  from  long-period 
surface  waves.  Street  et  al.  [1975]  empirically  determined  r0  =  100  km,  or  roughly 
twice  the  crustal  thickness.  Measuring  the  decay  rate  of  synthetic  Lg  phases  computed 
for  an  elastic  medium,  Herrmann  and  Kijko  [1983]  verified  that  Lg  frequency  domain 
spreading  was  accurately  described  as  cylindrical  and  substantiated  the  empirical  result 
of  Street  et  al.  [1975]  for  rQ.  Therefore  Lg  geometric  spreading  is  approximated  by 
(3.7)  with  r0  =  100  km  and  m  -  1/2. 


As  discussed  at  length  in  Section  1  of  this  report,  the  geometric  spreading  of  Pn  is  a 
complicated  function  of  the  velocity  gradient  in  the  upper  mantle,  and  it  is  probably 
frequency  dependent.  Therefore  a  simple  parameterization  such  as  (3.7)  is  not  likely 
to  be  applicable  to  Pn  over  a  broad  frequency  and  distance  range.  However,  over  lim¬ 
ited  distance  ranges  and  at  long  periods  a  power  law  range  dependence  may  be  ade¬ 
quate  to  approximate  the  spreading  of  Pn.  For  example.  Figure  3.2  plots  synthetic  Pn 
attenuation  at  1.25  Hz  for  an  elastic  earth  model  consisting  of  a  40-km-thick  homo¬ 
geneous  crust  over  a  mantle  with  a  compressional-wave  velocity  gradient  of  0.0013  s-1 
(Model  2,  Section  1).  This  gradient  corresponds  to  the  earth-flattening  transformation 
of  a  homogeneous  upper  mantle.  The  dashed  curve  is  the  best-fitting  power  law  model 
determined  by  least  squares.  The  spreading  rate  for  this  model  is  r-1 25.  Note  that  the 
difference  between  these  two  curves  is  considerably  less  than  the  stanJard  deviation  of 
typical  log  Af0  versus  magnitude  relations  (=  0.3)  [e.g.,  Veith  and  Clawson ,  1972; 
Dwyer  et  al.,  1983;  Hasegawa,  1983;  Sereno  et  al.,  1988].  Sereno  et  al.  [1988] 


33 


Log  Amplitude  (arbitrary  units) 


Pn  Geometric  Spreading 


Epicentral  Distance  (km) 


Figure  3.2.  Pn  geometric  spreading.  The  solid  curve  is  the  Pn  attenuation  with  dis¬ 
tance  at  1.25  Hz  calculated  for  an  elastic  earth  model  with  upper  mantle  velocity  gra¬ 
dient  equivalent  to  the  earth-flattening  transformation  of  a  homogeneous  medium 
(Model  2,  Section  2).  The  dashed  curve  plots  a  power  law  distance  dependence  r  23 
determined  by  a  least  squares  fit  to  the  synthetic  Pn  attenuation. 


estimate  Pn  Qlf)  from  data  recorded  at  NORESS  by  assuming  frequency-independent 
geometric  spreading  functions  with  range  dependencies  of  r~l,  r"lJ,  and  r-1-3.  For  each 
of  these  they  found  Qlf)  models  consistent  with  the  observed  spectra.  However,  r-13 
was  found  to  give  moments  that  were  consistent  with  the  moments  derived  from  the 
separate  Lg  inversion.  Thus,  the  Pn  spreading  model  near  1  Hz  derived  by  Sereno  et 
al.  [1988]  is  generally  consistent  with  an  upper  mantle  model  with  a  slight  positive 
velocity  gradient.  However,  the  Pn  Qlf)  determined  by  Sereno  et  al.  [1988]  includes 
any  frequency  dependence  occurring  in  the  geometric  spreading  function,  and  so  it 
probably  does  not  represent  upper  mantle  anelasticity. 

3.4.1  Lg  Inversion  Results 

The  Lg  inversion  includes  37  spectra  from  21  events  recorded  by  up  to  three  stations. 
The  spectra  were  inverted  between  0.5  and  2.5  Hz.  This  frequency  range  is  based  on 
a  comparison  of  the  Lg  spectra  to  the  pre-Lg  noise  spectra  normalized  to  the  signal 
window  length  (see  Appendix  B).  The  result  of  the  inversion  is  a  range  of  attenuation 
and  source  models  corresponding  to  a  broad  minimum  in  the  data  residuals.  For 
example.  Figure  3.3  plots  the  value  of  Lg  Q  at  1  Hz  (Qq)  against  the  exponent  of  a 
power  law  frequency  dependence  (q)  for  models  that  produce  data  variances  that  differ 
by  less  than  1%.  This  range  of  models  was  determined  by  fixing  q  and  inverting  for 
Go.  So  for  each  event,  and  the  comer  frequency  scaling  parameter.  In  general,  models 
with  low  Go  ^  associated  with  higher  source  levels  than  models  with  high  Go  since 
the  combination  these  parameters  is  constrained  to  fit  the  same  1-Hz  amplitude. 
Models  with  low  Go  are  also  associated  with  high  values  of  q  since  these  parameters 
are  constrained  by  spectral  shape.  Since  most  of  the  events  have  comer  frequencies 
greater  than  the  upper  limit  of  the  frequency  band  inverted,  the  attenuation  results  arc 
insensitive  to  the  source  assumptions. 

A  simple  constraint  that  can  be  applied  to  reduce  the  trade-off  among  model  parame¬ 
ters  is  to  require  consistency  in  the  source  levels  derived  for  each  event  from  data 
recorded  at  different  stations.  For  example,  Figure  3.4  plots  Lg  S0  for  two  of  the  Qlf) 
models  in  Figure  3.3.  The  model  with  higher  Q0  produces  more  consistent  source  lev¬ 
els  from  data  recorded  at  different  stations  than  the  model  with  lower  Go-  For  exam¬ 
ple,  the  source  levels  for  Events  9  and  17  derived  from  three  stations  are  nearly  the 
same  for  Lg  Qlf)  =  650  but  differ  by  over  50%  for  Lg  Qlf)  =  345/a\  The  difference  in 
epicentral  distance  between  the  furthest  and  closest  station  for  these  events  is  over  300 
km.  Events  1,  5,  11,  and  14  also  have  more  consistent  source  levels  for  the  higher  Go 
model.  In  general,  there  is  a  distance  dependence  in  the  derived  source  levels  for 
events  recorded  at  multiple  stations  for  models  with  Go  less  than  about  500.  This 
implies  that  Lg  Go  is  greater  than  or  equal  to  500.  Appendix  B  plots  the  Lg  spectra 
and  theoretical  spectra  based  on  the  model  with  Lg  Qlf)  »  650. 

Given  et  al.  [1989]  estimated  Lg  Q  by  applying  a  spectral  ratio  method  to  data 
recorded  from  one  of  the  10-ton  calibration  shots  (Event  20).  They  concluded  that  the 
average  value  of  Q  at  1  Hz  is  500  and  increases  with  frequency  although  the  frequency 
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Figure  3.3.  Lg  Q(f)  models  consistent  with  the  NRDC  observations.  The  value  of  Lg 
q  at  1  Hz  (Co)  is  plotted  against  the  exponent  of  a  power  law  frequency  dependence 
(q).  These  models  combine  with  different  source  levels  to  produce  data  variances  that 
differ  by  less  than  1%. 
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Figure  3.4.  Derived  source  levels  (S0)  for  (a)  Lg  Q(f)  =  650  and  (b)  Lg  Q(J)  =  345/ °4. 
The  event  numbers  are  from  Table  3.2.  Events  recorded  at  multiple  stations  have  an 
So  estimate  from  each  station.  Note  that  the  higher  Q0  model  reduces  the  variance  in 
the  derived  source  levels  from  multiple  station  observations. 


exponent  is  poorly  constrained.  They  approximate  the  frequency  dependence  as/04  to 
about  5  Hz  with  a  possible  increased  dependency  at  higher  frequency.  Their  Lg  Q  at  1 
Hz  is  consistent  with  the  value  obtained  here,  but  the  frequency  exponent  is  somewhat 
higher.  For  example,  the  model  in  Figure  3.3  with  Qq  equal  to  500  has  a  frequency 
exponent  close  to  0.2.  This  model  is  constrained  by  data  at  frequencies  less  than  2.5 
Hz.  A  comparison  of  the  pr e-Lg  noise  spectrum  with  the  Lg  signal  spectrum  for  this 
event  suggests  that  Sn  coda  may  contaminate  the  Lg  spectrum  at  KSU  beyond  4-5  Hz 
(Figure  A2.20).  Thus  the  higher  frequency  exponent  for  Lg  Q  determined  by  Given  et 
al.  [1989]  may  reflect  an  increased  contribution  of  Sn  coda  at  high  frequencies. 

The  earthquake  moments  are  estimated  from  Lg  S0  using  (3.5)  with  pc  =  2.5  gm/cm3 
and  pc  =  3.5  km/s.  The  S0  are  approximately  equal  to  the  logarithmic  mean  of  the 
source  level  below  the  comer  frequency.  Herrmann  and  Kijko  [1983]  use  synthetic  Lg 
spectra  to  demonstrate  that  the  logarithmic  mean  plus  one  standard  deviation  (close  to 
the  peak  rather  than  the  mean  spectral  amplitude)  should  be  used  for  the  S$  in  (3.5). 
Since  the  standard  deviation  of  log  S0  is  approximately  0.3,  the  S0  are  multiplied  by 
two  before  estimating  earthquake  moment  using  (3.5).  Figure  3.5  plots  earthquake 
moment  versus  mb  for  the  Lg  model  with  Q(f)  =  650.  The  symbols  indicate  the 
moments  of  the  five  events  in  Table  3.2  with  PDE  magnitudes  and  the  solid  lines  are 
the  moment-magnitude  relations  listed  in  Table  3.3.  Note  that  the  moments  derived  in 
this  study  are  generally  consistent  with  the  moment  versus  magnitude  relations  deter¬ 
mined  for  eastern  North  America  and  Scandinavia.  The  dashed  lines  in  Figure  3.5 
indicate  the  moments  of  earthquakes  that  are  not  listed  in  the  PDE  Bulletin  (Events  7, 
10,  and  17).  Comparing  these  to  the  moment  versus  magnitude  relations  indicates  that 
the  magnitudes  of  Events  7  and  17  are  between  4.1  and  4.4  and  that  the  magnitude  of 
Event  10  is  between  3.2  and  3.4. 

The  relationship  between  comer  frequency  (fc )  and  moment  for  the  Lg  Q  =  650  model 
is 


fc  =  6.26  x  107  \T0xr3  (3.8) 

which  implies  a  corner  frequency  of  2.6  Hz  for  the  largest  event  in  the  data  set  (Event 
5,  Table  3.2).  Since  this  is  greater  than  the  maximum  frequency  used  in  the  Lg  inver¬ 
sion,  the  comer  frequency  scaling  parameter  cannot  be  resolved  with  much  confidence. 
However,  Somerville  et  al.  [1987]  determined  a  similar  relation  for  events  in  eastern 
North  America.  They  determined  source  moment  and  duration  of  13  events  with  mNg 
between  4.7  and  6.6  by  modeling  waveforms  of  seismic  body  phases.  They  found  that 
the  relationship  between  source  duration  and  moment  is  consistent  with  constant  stress 
drop  scaling  and  that  the  median  stress  drop  is  approximately  100  bars.  Their  relation¬ 
ship  between  comer  frequency  and  moment  is 

fc  -  8.58  x  107  M-0m  (3.9) 

Note  that  both  (3.9)  and  (3.8)  give  comer  frequencies  greater  than  2.5  Hz  for  all  of  the 
events  studied. 
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Figure  3.5.  Moment  versus  magnitude  for  presumed  earthquakes  in  the  NRDC  data 
set.  The  moments  are  calculated  from  the  Lg  source  levels  using  (3.5)  and  are  plotted 
against  the  mb  reported  in  the  PDE  Bulletin  (Events  1,  2,  3,  5,  and  9  in  Table  3.2). 
The  solid  lines  are  moment  versus  magnitude  relations  determined  in  previous  studies 
(Table  3.3).  The  dashed  horizontal  lines  are  moments  of  earthquakes  not  reported  in 
the  PDE  Bulletin. 


Table  3.3  Moment  Versus  Magnitude  Relations 


Ml 

log  M0 

Region 

Reference 

-0.4-2.2 

0.90 Ml  +  17.50 

Scandinavia 

B ungum  et  al.  [1982] 

1. 8-4.9 

1.26AfL  +  16.44 

Central  United  States 

Dwyer  et  al.  [1983] 

1. 0-4.2 

0.94Aft  +  17.32 

E.  North  America 

Hasegawa  [1983] 

2. 1-4.3 

1.18A4  +  16.60t 

E.  North  America 

Shin  and  Herrmann  [1987] 

1. 1-3.8 

1.03A/L  +  17.10 

Scandinavia 

Sereno  et  al.  [1988] 

t  The  regression  coefficients  were  estimated  from  tabulated  moments  and  magnitudes. 


3.4.2  Pn  Inversion  Results 


The  Pn  inversion  includes  26  spectra  from  21  events  recorded  by  up  to  three  stations. 
The  spectra  were  inverted  between  1  and  10  Hz  based  on  signal-to-noise  ratio.  As 
was  the  case  for  Lg,  a  range  of  attenuation  and  source  models  produce  similar  data 
variances.  Figure  3.6  plots  Pn  Q0  against  r|  for  models  that  produce  data  variances  that 
differ  by  less  than  1%.  This  range  of  models  was  determined  by  fixing  r\  and  invert¬ 
ing  for  S0  for  each  event,  Q0,  and  the  comer  frequency  scaling  parameter. 

The  wide  range  of  Pn  Q(f)  models  in  Figure  3.6  cannot  be  distinguished  on  the  basis  of 
source  levels  derived  from  data  at  different  stations.  For  example,  Figure  3.7  plots  Pn 
SQ  for  two  very  different  Q(f)  models  that  fit  the  data  equally  well.  The  source  levels 
inferred  from  different  stations  are  no  more  consistent  for  one  of  these  models  than  for 
the  other.  Part  of  the  problem  is  that  there  are  only  a  few  events  with  Pn  data  from 
multiple  stations.  Other  factors  that  may  contribute  to  the  difficulty  are  (1)  an  invalid 
geometric  spreading  assumption,  (2)  radiation  pattern,  (3)  local  site  effects  and  (4)  pro¬ 
pagation  path  differences.  Whatever  the  cause,  the  consistency  of  derived  source  lev¬ 
els  cannot  be  used  to  resolve  trade-offs  among  source  and  path  parameters  as  was  done 
for  Lg.  For  this  reason.  Appendix  B  plots  the  Pn  spectra  and  theoretical  spectra  based 
on  derived  attenuation  models  that  represent  the  extremes  {Q(f)  =  1175  and  Q(f)  = 
300/ °5). 


The  explosion  moments  are  calculated  from  Pn  S0  using  (3.6)  with  ps  =  2.5  gm/cm3 
and  oq  =  5.0  km/s.  Events  13,  19,  and  20  are  known  to  be  explosions  [Thurber  et  al., 
1989].  Event  21  is  also  probably  an  explosion  in  the  Lake  Balkash  area  (see  Section 
3.2).  Other  likely  explosions  are  Events  11,  14,  16,  and  18,  based  on  low  LgIPn  ratios, 
origin  time  (morning  to  mid -afternoon),  location  and  general  aseismicity  of  the  region 
near  the  eastern  Kazakhstan  explosion  test  site  [Thurber  et  al.,  1989],  and  the  presence 
of  Rg.  The  Lg  magnitudes  of  these  events  are  calculated  using  the  correction  table 
applied  at  NORESS  (it  will  be  shown  in  Section  3.6  that  Lg  attenuation  in  eastern 
Kazakhstan  is  similar  to  that  observed  at  NORESS,  so  application  of  the  NORESS 
correction  table  is  reasonable).  Figure  3.8  plots  explosion  moment  versus  Lg  magni¬ 
tude  for  the  two  Q(f)  models.  The  moments  are  averaged  for  events  recorded  at  multi¬ 
ple  stations.  The  solid  line  is  the  moment-magnitude  relation  determined  by  Sereno  et 
al.  [1988]  for  mining  explosions  recorded  at  NORESS: 

log  Ml*p  =  1.08  Ml+  17.6  (3.10) 

The  model  with  Pn  Qif)  =  1175  gives  explosion  moments  consistent  with  (3.10),  while 
the  Q{f)  =  300f°5  model  fits  the  line  about  as  well  when  the  intercept  is  about  0.3 
higher.  Without  independent  estimates  of  source  moment  it  is  not  possible  to  deter¬ 
mine  which  of  these  models  more  accurately  separates  source  and  path  contributions  to 
the  observed  amplitudes. 
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Figure  3.6.  Pn  Q(f 
Q  at  1  Hz  (Co)  is 
Ol).  These  models 
differ  by  less  than 
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Pn  Q(f)  =  300f°-5 


Figure  3.7.  Derived  source  levels  (S0)  for  (a)  Pn  Q(f)  =  1175  and  ( b )  Pn  Q(f)  - 
30Qf05.  The  event  numbers  are  from  Table  3.2.  Events  recorded  at  multiple  stations 
have  an  So  estimate  from  each  station. 


Explosion  Moment  Versus  ML 


Figure  3.8.  Moment  versus  magnitude  for  presumed  explosions  in  the  NRDC  data  set. 
The  moments  are  calculated  from  the  Pn  source  levels  using  (3.6)  and  are  plotted 
against  Lg  magnitude  (Events  11,  13,  14,  16,  18,  19-21  in  Table  3.2).  The  asterisks 
are  moments  calculated  for  the  Pn  Q(f)  =  300f 05  model  and  the  squares  are  for  the  Pn 
Q{f)  =1175  model.  Note  that  the  moments  for  the  lower  Q0  model  are  about  a  factor 
of  2  higher  than  the  moments  associated  the  higher  Q0  model.  The  solid  line  is  the 
explosion  moment  versus  magnitude  relation  determined  by  Sereno  et  al.  [1988]  from 
NORESS  observations. 
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The  relationship  between  comer  frequency  and  explosion  moment  for  the  model  with 
PnQ(f)  =  1 175  is 


fe  =  4.65  x  107  Ml™  (3.11) 

This  relationship  for  the  Q(f)  =  3 OCf  0  5  model  is 

fc  =  6.31  x  107  Mlxri  (3.12) 

These  give  similar  comer  frequencies  since  the  M0  are  a  factor  of  two  higher  for  the 
lower  Q0  model.  Both  give  the  comer  frequency  of  an  ML  3.0  explosion  to  be  between 
5.0  and  6.0  Hz.  Note  that  Given  et  al.  [1989]  used  a  curve-fitting  algorithm  to  deter¬ 
mine  the  comer  frequency  of  Event  20  {ML  =  2.9)  to  be  5-6  Hz  from  Lg  spectra  and  7 
Hz  from  Pg  spectra. 

3.5  Noise  Spectra 

It  is  shown  in  this  section  that  ambient  noise  spectra  recorded  at  NORESS  are  similar 
to  noise  spectra  recorded  by  the  borehole  seismometer  at  the  NRDC  stations.  Berger 
et  al.  [1988]  provide  a  detailed  description  of  the  noise  at  the  NRDC  stations,  includ¬ 
ing  plots  of  ambient  noise  spectra  from  1  to  100  Hz  for  windy  and  calm  conditions  for 
the  borehole  and  surface  instruments  at  each  of  the  three  NRDC  stations.  This  section 
provides  a  direct  comparison  of  the  NORESS  noise  spectra  and  the  NRDC  station 
noise  spectra  computed  in  exactly  the  same  way  between  1  and  20  Hz. 

The  NORESS  single-channel  noise  spectrum  is  calculated  from  43  samples  with  5-s 
windows  [Henson  and  Bache,  1988].  Each  of  these  "samples"  is  the  average  spectrum 
for  all  elements  of  the  NORESS  array.  The  NRDC  station  noise  spectra  are  estimated 
from  samples  taken  prior  to  Pn  detections.  They  include  26  samples  with  5-s  win¬ 
dows,  and  data  from  the  three  NRDC  stations  have  been  averaged  together.  Figure  3.9 
plots  the  average  NRDC  ambient  noise  spectrum  computed  (solid  curve)  bounded  by 
one  standard  deviation  (dashed  curves).  The  dotted  curve  is  the  average  ambient  noise 
spectrum  at  NORESS  taken  from  Henson  and  Bache  [1988].  The  NORESS  noise 
spectrum  is  within  one  standard  deviation  of  the  average  noise  spectrum  for  the  NRDC 
stations. 

3.6  Comparison  to  NORESS  Data 

This  section  compares  the  parameterization  of  Lg  and  Pn  signal  spectra  derived  in  Sec¬ 
tion  3.4  for  the  NRDC  stations  to  the  parameterization  of  these  spectra  for  events 
recorded  at  NORESS  determined  by  Sereno  et  al.  [1988].  Figure  3.10  plots  theoretical 
Lg  spectra  at  800  km  for  an  ML  3.0  event  based  on  the  NORESS  and  NRDC  Q(f) 
models.  The  solid  curve  is  the  Lg  spectrum  for  one  of  the  models  ( Q(f)  =  650)  that  fits 
the  NRDC  station  data,  and  the  dashed  curve  is  the  Lg  spectrum  for  a  model  ( Q(f)  = 
560/ 026)  that  fits  the  NORESS  data.  The  same  source  spectrum  was  used  for  both 
curves.  The  spectra  are  very  similar  between  0.5  and  2.5  Hz  (the  frequency  band  used 
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Figure  3.9.  Ambient  noise  displacement  spectra  for  a  5-s  window.  The  solid  curve  is 
the  average  noise  spectrum  at  the  three  NRDC  stations  calculated  from  26  samples 
taken  prior  to  Pn.  The  dashed  curves  indicate  one  standard  deviation.  The  dotted 
curve  is  the  average  ambient  noise  spectrum  at  NORESS  based  on  43  samples  taken 
prior  to  Pn  [Henson  and  Bache,  1988]. 
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in  the  inversion  of  the  NRDC  data)  but  differ  at  higher  frequencies.  The  NORESS  | 

model  was  estimated  from  1  to  7  Hz  using  events  at  distances  greater  than  1000  km.  I 

At  these  ranges  the  higher-frequency  Sn  coda  influences  the  spectra  above  3  Hz,  and  . 

thus  raises  the  frequency  exponent  in  the  Q  model.  For  example,  Sereno  and  Bratt  1 

[1988]  show  that  NORESS  Lg  Q(f)  =  582/018  for  frequencies  between  1  and  3  Hz. 

Thus,  the  difference  between  the  two  Lg  spectra  above  3  Hz  in  Figure  3.10  may  be  < 

due  to  the  effect  of  Sn  coda  included  in  the  NORESS  spectra,  but  not  seen  in  the  lower  I 

frequency  NRDC  station  spectra. 

A  range  of  Pn  attenuation  and  source  models  was  determined  in  Section  3.4,  and  it  is  | 

not  possible  to  resolve  the  parameter  trade  offs  without  independent  estimates  of 
source  moment.  Therefore,  in  Figure  3.11  Pn  spectra  for  two  extreme  models  con-  I 

sistent  with  the  NRDC  spectra  from  1  to  10  Hz  are  compared  to  Pn  spectra  for  the  I 

NORESS  model  determined  by  Sereno  et  al.  [1988].  These  model  Pn  spectra  are  com¬ 
puted  at  800  km  for  an  ML  3.0  explosion.  Table  3.4  lists  the  parameters  for  the  calcu-  I 

lation.  The  two  solid  curves  representing  the  NRDC  data  are  very  similar  for  the  1-10  I 

Hz  frequency  band  used  in  the  inversion.  The  NORESS  spectrum  has  a  similar  shape 
between  1  and  10  Hz,  but  its  amplitude  is  lower  by  about  a  factor  of  two.  The  I 

differences  beyond  10  Hz  are  caused  by  differences  in  the  source  parameters.  The  I 

comer  frequency  of  the  NORESS  model  for  an  ML  3.0  explosion  is  10  Hz,  while  for 
the  NRDC  models  it  is  between  5  and  6  Hz.  However,  the  attenuation  model  at  j 

NORESS  is  included  in  the  range  of  attenuation  models  found  to  fit  the  NRDC  data  J 

between  1  and  10  Hz  (compare  the  first  and  last  rows  in  Table  3.4).  The  main 
difference  is  that  the  Pn  amplitude  is  smaller  at  NORESS  for  fixed  source  moment. 

3.7  Summary 

Regional  Pn  and  Lg  spectra  recorded  by  three  seismic  stations  near  the  eastern  Kaza¬ 
khstan  nuclear  explosion  test  site  are  inverted  for  source  and  attenuation  parameters. 

These  results  am  "•'mpared  to  results  obtained  with  the  same  methods  using  Pn  and  Lg 
spectra  recorded  by  the  NORESS  array  in  Scandinavia.  The  main  conclusions  are: 

1.  The  attenuation  of  Lg  spectra  between  0.5  and  2.5  Hz  in  eastern  Kazakhstan  is  j 

consistent  with  cylindrical  geometric  spreading  and  a  constant  Q  equal  to  650. 
Frequency-dependent  models  are  also  consistent  with  the  Lg  spectra.  However,  these 
models  produce  a  greater  variance  in  the  moments  estimated  for  the  same  event  J 

recorded  at  multiple  stations. 

2.  The  Pn  spectra  between  1  and  10  Hz  in  eastern  Kazakhstan  are  consistent  with 
r  ~13  geometric  spreading  and  a  wide  range  of  Q(f)  models  varying  from  a  constant  Q 
of  1175  to  a  Q<S)  =  22Q f0b:  Each  of  these  models  provides  a  different  set  of  source 
moment  estimates,  but  an  indepenent  estimate  of  source  moment  is  necessary  to  favor 
one  and  resolve  the  parameter  trade  offs.  It  should  be  emphasized  that  the  Pn  Q 
models  include  the  frequency  dependence  of  geometric  spreading,  and  therefore  they 
are  not  a  direct  estimate  for  upper  mantle  anelasticity. 
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Figure  3.11.  Theoretical  Pn  spectra  at  800  km  for  a  magnitude  3.0  explosion  based  on 
the  inversion  results  using  the  NRDC  (solid  curves)  and  NORESS  data  (dashed).  The 
results  for  two  NRDC  Pn  models  are  plotted  since  parameter  trade-offs  cannot  be 
resolved  on  the  basis  of  source  moment.  Model  parameters  are  listed  in  Table  3.4. 
The  bottom  curve  is  the  average  NORESS  ambient  noise  scaled  to  a  10-s  window 
length. 
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Table  3.4  Pn  Attenuation  and  Source  Parameters 


Database 

Qo 

n 

log  M0 

C=fc<ri 

NORESS 

325 

0.48 

1.08  Ml  +  17.6 

8.77  xlO7 

NRDC 

1175 

0.00 

1.08A/,,  +  17.6 

4.65  xlO7 

NRDC 

300 

0.50 

1.08A/l  +  17.9 

6.31  xlO7 
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3.  The  average  ambient  noise  spectrum  at  NORESS  falls  within  one  standard  devia¬ 
tion  of  the  average  noise  spectrum  observed  at  the  three  NRDC  stations  prior  to  the  26 
Pn  recordings  used  in  this  study. 

4.  Lg  attenuation  near  the  nuclear  explosion  test  site  in  eastern  Kazakhstan  is  very 
similar  to  the  Lg  attenuation  observed  at  NORESS  in  the  0.5  to  2.5  Hz  band.  Extrapo¬ 
lation  of  the  NRDC  station  model  to  higher  frequency  indicates  differences  in  propaga¬ 
tion  characteristics  but  this  difference  may  simply  reflect  different  contributions  from 
Sn  coda  because  of  the  different  frequency  bands  used  to  determine  the  different  Q 
models. 

5.  The  Pn  spectral  shape  in  eastern  Kazakhstan  is  similar  to  that  observed  at 
NORESS  between  1  and  10  Hz.  The  attenuation  model  derived  by  Sereno  et  al. 
[1988]  for  NORESS  (£)(/)  =  325/ °'48)  fits  the  NRDC  station  data  as  well  as  any  of  the 
other  models.  However,  the  the  Pn  amplitudes  between  1  and  5  Hz  are  about  a  factor 
of  two  higher  at  the  NRDC  stations  than  at  NORESS  for  the  same  source  moment. 
Since  the  noise  levels  in  the  two  regions  are  nearly  the  same,  this  implies  that  the 
single-station  magnitude  thresholds  for  detecting  Pn  appear  to  be  as  much  as  0.3  lower 
in  eastern  Kazakhstan  than  they  are  at  NORESS. 
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4.  CONCLUSIONS 


The  overall  goal  of  this  work  is  to  estimate  the  detection  and  location  capability  of 
actual  and  hypothetical  networks  of  stations  and  arrays  in  and  around  the  Soviet 
Union.  These  estimates  require  the  simulation  of  the  detection  and  location  capability 
of  the  networks.  This  report  summarizes  two  studies  of  important  issues  involved  in 
the  normalization  of  the  simulations.  The  first  is  a  theoretical  study  of  the  dependence 
of  Pn  geometric  spreading  on  the  velocity  gradient  in  the  upper  mantle.  This  is  impor¬ 
tant  because  most  methods  for  estimating  Pn  attenuation  require  assumptions  regarding 
geometric  spreading.  The  second  study  is  an  empirical  parameterization  of  regional 
wave  spectra  recorded  near  the  eastern  Kazakhstan  nuclear  explosion  test  site.  The 
results  arc  compared  to  results  from  a  similar  parameterization  of  data  recorded  by  the 
NORESS  array  in  Norway.  Data  recorded  by  the  NORESS  and  ARCESS  arrays  are 
plentiful  and  provide  an  excellent  basis  for  normalizing  the  simulations.  Data  recorded 
in  eastern  Kazakhstan  are  relatively  few  at  this  time,  but  can  be  used  to  see  if  the 


NORESS/ARCESS  data  can  be  used  as  a  reasonable  representation  of  conditions  in. 


this  area 
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The  main  Conclusions  of  the  first  study  are: 
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•  Pn  geometric  spreading  is  extremely  sensitive  to  the  velocity  gradient  in  the  upper 
mantle.  For  example,  the  Pn  amplitude  at  1000  km  for  a  model  with  the  weak 
velocity  gradient  corresponding  to  the  earth-flattening  transformation  of  a  homo¬ 
geneous  upper  mantle  is  a  factor  of  20  larger  than  the  Pn  amplitude  for  a  layer 
over  a  halfspace. 


•  Pn  geometric  spreading  for  an  upper  mantle  model  with  a  linear  velocity  gradient 
is  frequency-dependent.  The  Pn  waveform  for  a  model  with  a  weak  velocity  gra¬ 
dient  is  similar  to  an  head  wave  at  distances  within  a  few  hundred  kilometers  of 
the  critical  distance,  but  has  the  time  dependence  of  a  turning  wave  at  longer 
ranges  (the  derivative  of  the  head  wave). 


•  Methods  that  assume  frequency-independent  Pn  geometric  spreading  (e.g.,  spectral 
ratio  methods)  can  lead  to  overestimates  of  upper  mantle  Q  if  the  velocity  in  the 
upper  mantle  increases  with  depth. 

These  results  have  important  implications  for  the  interpretation  of  Pn  amplitudes  in 
terms  of  elastic  and  anelastic  upper  mantle  structure  and  for  practical  issues  related  to 
regional  seismic  monitoring  of  nuclear  explosion  testing. 

The  parameterization  of  regional  wave  spectra  used  in  the  second  study  requires 
assumption  of  the  geometric  spreading.  The  Pn  spectra  recorded  in  eastern  Kazakhstan 
are  parameterized  by  a  frequency-independent  geometric  spreading  function  with  r  “u 
distance  dependence.  This  is  the  approximate  distance  dependence  at  1  Hz  between 
300  and  and  1000  km  for  a  homogeneous  upper  mantle  (i.e.,  the  gradient  of  the  earth- 
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flattening  transformation).  Since  the  spreading  is  assumed  to  be  independent  of  fre¬ 
quency,  the  Pn  Q(f)  models  include  frequency-dependent  effects  associated  with  the 
velocity  structure  of  the  upper  mantle  and  do  not  have  a  direct  interpretation  in  terms 
of  upper  mantle  anelasticity.  However,  a  distinction  between  anelastic  and  geometric 
effects  is  not  required  to  represent  the  distance  dependence  of  spectra. 

The  main  conclusions  of  the  second  study  are: 

•  The  attenuation  of  Lg  spectra  between  0.5  and  2.5  Hz  in  eastern  Kazakhstan  is 
nearly  identical  to  the  Lg  attenuation  in  Scandinavia. 

•  The  Pn  spectra  between  1  and  10  Hz  in  eastern  Kazakhstan  are  similar  in  shape 
to  the  Pn  spectra  in  Scandinavia  at  similar  distances.  However,  the  absolute  lev¬ 
els  appear  to  be  as  much  as  a  factor  of  two  higher  in  eastern  Kazakhstan  for  fixed 
source  moment  (given  that  source  moment  is  related  to  Lg  magnitude). 

•  The  average  ambient  noise  spectrum  at  NORESS  falls  within  one  standard  devia¬ 
tion  of  the  average  noise  spectrum  at  the  three  NRDC  stations. 

It  follows  from  the  above  conclusions  that  the  magnitude  thresholds  for  detecting  Pn 
appear  to  be  as  much  as  0.3  higher  in  Scandinavia  than  in  eastern  Kazakhstan.  Since 
Lg  propagation  is  nearly  the  same,  the  magnitude  thresholds  for  detecting  either  Pn  or 
Lg  are  less  than  0.3  higher  in  Scandinavia  than  in  eastern  Kazakhstan.  It  is  uncertain 
which  of  these  Pn  models  more  accurately  describes  the  average  attenuation  in  the 
Soviet  Union,  but  the  comparison  gives  a  measure  of  the  amount  of  regional  variability 
that  may  be  expected. 
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APPENDIX  A:  NRDC  WAVEFORMS 


This  appendix  plots  the  waveform  data  for  all  events  used  in  the  regional  wave 
attenuation  study  in  eastern  Kazakhstan  described  in  Section  3  of  this  report.  The  sta¬ 
tion  identification  and  epicentral  distance  are  listed  in  the  upper  right  comer  of  each 
plot.  Only  vertical-component  recordings  from  the  borehole  sensor  at  each  station  are 
plotted.  The  seismograms  are  filtered  between  0.5  and  20  Hz.  The  amplitudes  are 
converted  to  nanometers  from  digital  counts  using  calibration  values  at  1  Hz  obtained 
from  the  Center  for  Seismic  Studies  (these  calibration  values  are  listed  in  each  figure 
caption).  Note  that  these  are  only  approximate  amplitudes  since  the  instrument 
response  (flat  to  velocity  between  0.2  and  100  Hz)  has  not  been  deconvolved.  Pn  and 
Lg  phases  are  marked  in  each  figure  with  a  vertical  line. 
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2.  Waveform  for  Event  2  (Table  3.2).  This  is  an  mb  4.3  earthquake 
n  the  PDE  Bulletin.  The  calibration  value  is  0.012929  nanometers  per  digital 
station  KKL. 


Waveform  for  Event  3  (Table  3.2).  This  is  an  mb  4.6  earthquake 
PDE  Bulletin.  The  calibration  value  is  0.012929  nanometers  per  digital 


Figure  A. 4.  Waveform  for  Event  4  (Table  3.2).  The  event  type  is  unknown.  The 

calibration  value  is  0.012929  nanometers  per  digital  count  for  station  KSU. 


Figure  A.5.  Waveforms  for  Event  5  (Table  3.2).  This  is  an  mb  4.6  earthquake 
reported  in  the  PDE  Bulletin.  The  calibration  values  are  0.411528  and  1.638321 
nanometers  per  digital  count  for  stations  KKL  and  BAY,  respectively. 


Figure  A. 6.  Waveforms  for  Event  6  (Table  3.2).  The  event  type  is  unknown.  The 
calibration  values  are  0.012781  and  0.012929  nanometers  per  digital  count  for  stations 
KKL  and  BAY,  respectively. 


Waveform  for  Event  7  (Table  3.2).  This  is  probably  an  earthquake  but 
ted  in  the  PDE  Bulletin.  The  calibration  value  is  0.012929  nanometers 
nt  for  station  BAY. 


calibration  value  is  0.006406  nanometers  per  digital  count  for  station  KKL. 


Figure  A.9.  Waveforms  for  Event  9  (Table  3.2).  This  is  an  mb  4.4  earthquake 
reported  in  the  PDE  Bulletin.  Note  the  exceptionally  strong  Sn  at  KKL  and  BAY. 
The  calibration  value  is  0.006406  nanometers  per  digital  count  for  all  three  stations. 


Event  11 


ire  A.  11.  Waveforms  for  Event  11  (Table  3.2).  This  is  probably  an  explosion. 
Lg  magnitude  is  approximately  2.6.  Note  the  small  LgIPn  ratio  for  this  event, 
calibration  value  is  0.006406  nanometers  per  digital  count  for  both  stations. 


Figure  A.  12.  Waveform  for  Event  12  (Table  3.2).  The  event  type  is  unknown.  The 
calibration  value  is  0.006406  nanometers  per  digital  count  for  KSU. 


Event  14 
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Figure  A.  14.  Waveforms  for  Event  14  (Table  3.2).  This  is  probably  an  explosion. 
The  Lg  magnitude  of  this  event  is  approximately  2.3.  The  calibration  value  is 
0.006406  nanometers  per  digital  count  at  both  stations. 


m 


Figure  A.  16.  Waveform  for  Event  16  (Table  3.2).  The  event  is  probably  an  explo¬ 
sion.  The  Lg  magnitude  for  this  event  is  approximately  2.6.  The  calibration  value  is 
0.006406  nanometers  per  digital  count  for  station  KSU. 


Event  17 


Waveforms  for  Event  17  (Table  3.2).  This  event  is  an  earthquake  in 
irea  [Thurber  et  al.,  1989].  The  calibration  value  is  0.006406  nanome- 
;ount  for  all  three  stations. 


Figure  A.  18.  Waveforms  for  Event  18  (Table  3.2).  This  event  is  probably  an  explo¬ 
sion.  The  Lg  magnitude  is  approximately  2.5.  The  calibration  value  is  0.006406 
nanometers  per  digital  count  for  both  stations. 


Event  19 
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Figure  A.  19.  Waveforms  for  Event  19  (Table  3.2).  This  event  is  a  10-ton  calibration 
explosion  [Thurber  et  al„  1989;  Given  et  al.,  1989].  The  vertical  scale  is  adjusted  for 
the  waveform  at  KSU  (the  data  at  BAY  and  KKL  are  not  clipped).  The  Lg  magnitude 
for  this  event  is  approximately  2.9.  The  calibration  value  for  station  KSU  is  0.006406 
nanometers  per  digital  count  and  it  is  1.638321  for  stations  BAY  and  KKL. 


Event  21 
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Figure  A.21.  Waveforms  for  Event  21  (Table  3.2).  This  event  is  an  explosion  in  the 
Lake  Balkash  area  (compare  to  Event  13).  Note  the  strong  Rg  phase.  The  Lg  magni¬ 
tude  for  this  event  is  approxitrr  :ely  2.9.  The  calibration  value  is  0.012929  nanometers 
per  digital  count  for  station  BAY. 


APPENDIX  B:  NRDC  SPECTRA  AND  INVERSION  RESULTS 


I 


Regional  Pn  and  Lg  spectra  from  21  events  (Table  3.2)  were  inverted  for  Q(J),  long- 
period  source  level  (S0),  and  a  scaling  parameter  (c)  that  relates  comer  frequency  and 
S0  for  all  events  in  the  data  set.  Earthquake  moment  is  related  to  Lg  S0  by  (3.5)  and 
explosion  moment  is  related  to  Pn  S0  by  (3.6).  The  geometric  spreading  was  modeled 
as  (3.7)  with  r0  -  100  km  and  m  =  1/2  for  Lg  and  r0  =  1  km  and  m  =  1.3  for  Pn.  The 
path  result  for  Lg  between  0.5  and  2.5  Hz  is  Q(f)  =  650.  The  Pn  inversion  between  1 
and  10  Hz  produced  a  range  of  Q(f)  and  source  models  that  fit  the  observed  data  with 
nearly  equal  fidelity.  The  range  of  attenuation  models  is  approximately  Q(f)  =  1175  to 
Q(f)  =  300/ °5. 

Table  B.l  lists  the  source  parameters  for  the  model  with  Pn  Q(f)  =  1175.  The  mb  are 
from  the  PDE  Bulletin.  The  ML  were  calculated  from  the  Lg  amplitude  and  the 
NORESS  correction  table.  Table  B.2  lists  the  source  parameters  for  the  model  with  Pn 
Q(f)  =  300/ °5.  Table  B.3  lists  the  source  parameters  derived  from  the  Lg  inversion. 
Note  that  comer  frequency  is  not  a  parameter  of  the  inversion  for  each  event,  but  is 
calculated  from  the  long-period  level  and  the  scaling  parameter.  Since  the  Lg  inver¬ 
sion  only  included  frequencies  less  than  2.5  Hz,  comer  frequency  is  constrained  only 
by  a  small  number  of  large  events. 

Figures  B.1-B.21  display  theoretical  and  data  spectra  for  all  of  the  events  used  in  the 
inversion.  The  theoretical  spectra  are  computed  using  (3.4)  with  the  inversion  results 
for  Q(f)  and  S(f).  Two  theoretical  spectra  are  plotted  for  each  Pn  observation 
corresponding  to  the  models  in  Tables  B.l  and  B.2.  The  more  coarsely  dashed  of  the 
two  curves  corresponds  to  the  model  in  Table  B.2  ( Q(f)  =  300/  °'5).  The  lowest  curve 
in  each  Pn  plot  is  the  ambient  noise  spectrum  calculated  from  a  5-s  window  starting 
5.5  s  before  Pn.  The  theoretical  Lg  spectra  are  calculated  from  the  parameters  in  Table 
B.3.  The  lowest  curve  is  the  pre-Lg  noise  spectrum  calculated  from  a  5-s  window 
starting  5.5  s  before  Lg.  The  Pn  and  Lg  noise  spectra  are  scaled  to  the  signal  window 
length  using  (3.1).  The  Pn  spectra  are  plotted  from  0.5  to  20  Hz  and  the  Lg  spectra 
are  plotted  from  0.5  to  5.0  Hz.  The  spectra  are  ordered  for  each  event  with  data  from 
the  closest  station  plotted  at  the  top  of  each  figure.  The  Pn  spectra  are  plotted  on  the 
left  and  the  Lg  spectra  are  plotted  on  the  right. 
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Table  B.l  Pn  Inversion  Results  {Q(j)  =  1175) 


Event 

mb 

Ml 

Type 

Station 

Distance 

(km) 

So 

(cm2-s) 

fc 

(Hz) 

1 

4.3 

EQ 

2 

4.3 

EQ 

3 

4.6 

EQ 

4 

KSU 

642 

810.5 

3.2 

5 

4.6 

EQ 

KKL 

742 

391.5 

4.0 

6 

KKL 

378 

9.6 

13.9 

7 

EQ 

BAY 

977 

268.3 

4.6 

8 

KKL 

499 

41.0 

8.6 

9 

4.4 

EQ 

BAY 

1047 

335.3 

4.2 

KKL 

956 

516.3 

3.7 

KSU 

724 

162.2 

5.4 

10 

EQ 

BAY 

934 

296.4 

4.4 

11 

2.6 

EX 

BAY 

386 

103.2 

6.3 

KKL 

399 

47.8 

8.1 

12 

KSU 

640 

108.0 

6.2 

13 

2.5 

EX 

BAY 

456 

19.9 

10.9 

KKL 

309 

31.8 

9.3 

KSU 

436 

16.8 

11.5 

14 

2.3 

EX 

BAY 

555 

34.6 

9.1 

15 

KKL 

492 

107.0 

6.2 

KSU 

242 

10.9 

13.3 

16 

2.6 

EX 

KSU 

644 

59.8 

7.5 

17 

EQ 

KKL 

996 

1042.6 

2.9 

KSU 

777 

548.8 

3.6 

18 

2.5 

EX 

BAY 

785 

57.3 

7.7 

KSU 

577 

64.4 

7.4 

19 

2.9 

EX 

KSU 

636 

116.4 

6.0 

20 

2.9 

EX 

KSU 

636 

150.4 

5.5 

21 

2.9 

EX 

BAY 

477 

54.3 

7.8 
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Table  B.2  Pn  Inversion  Results  (Q(f)  =  300/ °5) 


Event 

mb 

Ml 

Type 

Station 

Distance 

(km) 

So 

(cm2-s) 

fc 

(Hz) 

1 

4.3 

EQ 

2 

4.3 

EQ 

3 

4.6 

EQ 

4 

KSU 

642 

1686.0 

3.4 

5 

4.6 

EQ 

KKL 

742 

947.5 

4.1 

6 

KKL 

378 

15.6 

16.0 

7 

EQ 

BAY 

977 

913.3 

4.1 

8 

KKL 

499 

74.6 

9.5 

9 

4.4 

EQ 

BAY 

1047 

1260.4 

3.7 

KKL 

956 

1755.1 

3.3 

KSU 

724 

388.5 

5.5 

10 

EQ 

BAY 

934 

956.1 

4.1 

11 

2.6 

EX 

BAY 

386 

161.8 

7.3 

KKL 

399 

77.8 

9.4 

12 

KSU 

640 

230.7 

6.5 

13 

2.5 

EX 

BAY 

456 

35.3 

12.2 

KKL 

309 

45.9 

11.2 

KSU 

436 

29.3 

13.0 

14 

2.3 

EX 

BAY 

555 

67.7 

9.8 

15 

KKL 

492 

193.8 

6.9 

KSU 

242 

15.2 

16.2 

16 

2.6 

EX 

KSU 

644 

126.3 

8.0 

17 

EQ 

KKL 

996 

3824.7 

2.6 

KSU 

777 

1426.8 

3.6 

18 

2.5 

EX 

BAY 

785 

146.3 

7.6 

KSU 

577 

127.2 

8.0 

19 

2.9 

EX 

KSU 

636 

253.4 

6.3 

20 

2.9 

EX 

KSU 

636 

321.7 

5.8 

21 

2.9 

EX 

BAY 

477 

96.7 

8.7 
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Table  B.3  Lg  Inversion  Results  (Q(J)  =  650) 


Event 

mb 

Ml 

Type 

Station 

Distance 

(km) 

So 

(cm2-s) 

fc 

(Hz) 

1 

4.3 

EQ 

BAY 

1283 

2711.2 

3.2 

KKL 

1118 

2530.5 

3.3 

2 

4.3 

EQ 

KKL 

1069 

2230.1 

3.4 

3 

4.6 

EQ 

KKL 

1087 

1809.6 

3.7 

4 

KSU 

642 

1205.6 

4.2 

5 

4.6 

EQ 

BAY 

898 

5254.0 

2.6 

KKL 

742 

4786.9 

2.7 

6 

BAY 

484 

54.4 

11.9 

KKL 

378 

48.0 

12.4 

7 

EQ 

BAY 

977 

1491.6 

3.9 

8 

KKL 

499 

122.4 

9.1 

9 

4.4 

EQ 

BAY 

1047 

1908.6 

3.6 

KKL 

956 

2239.6 

3.4 

KSU 

724 

1778.2 

3.7 

10 

EQ 

BAY 

934 

304.5 

6.7 

KKL 

825 

276.9 

6.9 

11 

2.6 

EX 

BAY 

386 

47.6 

12.4 

KKL 

399 

47.0 

12.5 

12 

13 

2.5 

EX 

BAY 

456 

54.4 

11.9 

KKL 

309 

41.2 

13.0 

KSU 

436 

32.3 

14.1 

14 

2.3 

EX 

BAY 

555 

29.2 

14.6 

KSU 

199 

28.0 

14.8 

15 

KKL 

492 

205.8 

7.6 

KSU 

242 

104.5 

9.6 

16 

2.6 

EX 

KSU 

644 

48.6 

12.3 

17 

EQ 

BAY 

1105 

1549.1 

3.9 

KKL 

996 

1575.5 

3.9 

KSU 

777 

1472.3 

4.0 

18 

2.5 

EX 

KSU 

577 

74.7 

10.7 

19 

2.9 

EX 

BAY 

246 

63.5 

11.3 

KKL 

253 

116.5 

9.2 

KSU 

636 

62.8 

11.3 

20 

2.9 

EX 

BAY 

246 

72.1 

10.8 

KKL 

253 

141.4 

8.6 

KSU 

636 

73.6 

10.7 

21 

2.9 

EX 

BAY 

477 

77.6 

10.6 
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Figure  B.4.  Spectra  for  Event  4  (Table  3.2).  See  text  for  explanation. 
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Figure  B.5.  Spectra  for  Event  5  (Table  3.2).  See  text  for  explanation. 
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Figure  B.6.  Spectra  for  Event  6  (Table  3.2).  See  text  for  explanation. 
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Figure  B.8.  Spectra  for  Event  8  (Table  3.2).  See  text  for  explanation. 
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Figure  B.16.  Spectra  for  Event  16  (Table  3.2). 
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Figure  B.17.  Spectra  for  Event  17  (Tabic  3.2).  See  text  for  explanation. 
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I  Figure  B.21.  Spectra  for  Event  21  (Table  3.2).  See  text  for  explanation. 
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